


‘ 


‘W S aN 
iby F ‘ hh Me , 









Institutional Archive of the Naval Postgraduate School 


Calhoun: The NPS Institutional Archive 
DSpace Repository 


Theses and Dissertations 1. Thesis and Dissertation Collection, all items 


1997-06 


Nearshore currents over a barred beach 


Faria, Garcez; Fernando, Antonio 


Monterey, California. Naval Postgraduate School 


http://ndl.handle.net/10945/8729 


Downloaded from NPS Archive: Calhoun 


| Calhoun is the Naval Postgraduate School's public access digital repository for 
(8 D U DLEY research materials and institutional publications created by the NPS community. 
«ist | | Calhoun is named for Professor of Mathematics Guy K. Calhoun, NPS's first 


lil \ KNOX appointed — and published -—- scholarly author. 


http://www.nps.edu/library 






LIBRARY Dudley Knox Library / Naval Postgraduate School 
411 Dyer Road / 1 University Circle 
Monterey, California USA 93943 


NPS ARCHIVE 
1997, Oc 
GARCEZ FARIA, A. 


NAVAL POSTGRADUATE SCHOOL 
Monterey, California 





DISSERTATION 


NEARSHORE CURRENTS OVER 
_A BARRED BEACH 


by 


Antonio Fernando pC Faria 
} 


(mmew1 997 
Dissertation Supervisor and 
Committee Chairman: Edward B. Thornton 
Committee Members: Timothy P. Stanton 
_— Thomas H.C. Herbers 
Thesis Roger T. Williams 
G18053 Thomas C. Lippmann 


Approved for public release; distribution is unlimited. 





DUDLEY KNOX LIBRARY 
NAW. POSTGRADUNTE. ScHOoL 
MOMMEREY CA Sofia 34 





REPORT DOCUMENTATION PAGE Form Approved OMB No. 0704-0188 | 


f 


"| 


| Public reporung burden for this collection of informauon 1s esumated to average | hour per response, including the time for reviewing instruction, searching existing data 

| sources, gathering and maintaining the data needed, and completing and reviewing the collection of information. Send comments regarding this burden estimate or any 

f other aspect of ths collection of information, including suggestions tor reducing Wis buiden, to Washiugtua Headquarters Services, Directorate for Information 

| Operations and Reports, 1215 Jefferson Davis Highway, Suite 1204, Arlington, VA 22202-4302, and to the Office of Management and Budget, Paperwork Reduction 
Project (0704-0188) Washington DC 20503. 


1. AGENCY USE ONLY (Leave blank) |2. REPORT DATE 3. REPORT TYPE AND DATES COVERED 
June 1997 Doctoral Dissertation 
4. TITLE AND SUBTITLE NEARSHORE CURRENTS OVER A BARRED FUNDING NUMBERS 
BEACH 


6. AUTHOR(S) Faria, Antonio Fernando Garcez 


7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) 
Naval Postgraduate School 
Monterey CA 93943-5000 


9. SPONSORING/MONITORING AGENCY NAME(S) AND ADDRESS(ES) 
Office of Naval Research 
Coastal Sciences 





/ 


| 











8. PERFORMING 
ORGANIZATION 
REPORT NUMBER 
























SPONSORING/MONITORING 
AGENCY REPORT NUMBER 







11. SUPPLEMENTARY NOTES The views expressed in this thesis are those of the author and do not reflect the 
official policy or position of the Department of Defense or the U.S. Government. 


12a. DISTRIBUTION/AVAILABILITY STATEMENT 12b. DISTRIBUTION CODE 
Approved for public release; distribution is unlimited. 


13. ABSTRACT (maximum 200 words) 

The objective of this dissertation is to develop numerical models and compare their predictions with data 
acquired during the DUCK94 experiment in order to improve our physical understanding of the hydrodynamic 
processes governing the vertical and cross-shore distnibutions of both longshore and cross-shore currents over < 
barred beach. The vertical structure of the mean longshore current is found to be well described by a logarithmic 
profile and a relationship between bed shear stress and bottom roughness, including the influence of mpples and 
mega-ripples, was also found. The vertical structure of the mean cross-shore current (undertow) is modeled using 
an eddy viscosity closure scheme to solve for the turbulent shear stress and includes contributions from breaking 
wave rollers (Lippmann and Thornton, 1997). These models of the vertical profiles of longshore and cross-shore 
mean currents are combined to formulate a quasi three-dimensional model to describe the cross-shore distribution 
of the longshore current. This model includes turbulent mixing due to the cross-shore advection of mean momentum 
of the longshore current by the mean cross-shore current (Putrevu and Svendsen, 1993) and contnibutions from wave 
rollers. 














14. SUBJECT TERMS longshore current, nearshore, undertow, radiation stress, bottom 15. NUMBER OF 


stress, momentum mixing. PAGES 178 
16. PRICE CODE 
)“=— ss; 

17. SECURITY CLASSIFI- 18. SECURITY CLASSIFI- 19. SECURITY CLASSIFICA- | 20. LIMITATION OF 
CATION OF REPORT CATION OF THIS PAGE TION OF ABSTRACT ABSTRACT 
Unclassified Unclassified Unclassified UL 

NSN 7540-01-280-5500 Standard Form 298 (Rev. 2-89) 


Prescribed by ANSI Std. 239-18 298-102 _ 


i] 





Approved for public release; distribution is unlimited. 


NEARSHORE CURRENTS OVER 
A BARRED BEACH 


Antonio Fernando Garcez Faria 
Lieutenant Commander, Brazilian Navy 
B.S., Brazilian Naval Academy, 1983 
M.S., Naval Postgraduate School, 1995 


Submitted in partial fulfillment of the requirements for the degree of 
DOCTOR OF PHILOSOPHY IN PHYSICAL OCEANOGRAPHY 
from the 


NAVAL POSTGRADUATE SCHOOL 


June eS ers 


\ Ps HEC W(Vs | 
(ACT, Oo Gy, 
OARCEZ FARIA CF 


DUDLEY KNOX LIBRARY 
TGRADUATE SCHOOL 


~~ 


= CG ‘ a3 wv 


ABSTRACT 


The objective of this dissertation is to develop numerical models and compare their 
predictions with data acquired during the DUCK94 experiment in order to improve our 
physical understanding of the hydrodynamic processes governing the vertical and cross-shore 
distributions of both longshore and cross-shore currents over a barred beach. The vertical 
structure of the mean longshore current is found to be well described by a logarithmic profile 
and a relationship between bed shear stress and bottom roughness, including the influence 
of ripples and mega-ripples, was also found. The vertical structure of the mean cross-shore 
current (undertow) is modeled using an eddy viscosity closure scheme to solve for the 
turbulent shear stress and includes contributions from breaking wave rollers (Lippmann and 
Thornton, 1997). These models of the vertical profiles of longshore and cross-shore mean 
currents are combined to formulate a quasi three-dimensional model to describe the cross- 
shore distribution of the longshore current. This model includes turbulent mixing due to the 
cross-shore advection of mean momentum of the longshore current by the mean cross-shore 


current (Putrevu and Svendsen, 1993) and contributions from wave rollers. 
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I. PREFACE 

The nearshore currents generated by obliquely incident breaking waves within the 
surf zone can exceed 1 ms” during high energy conditions, and are some of the strongest 
currents with the greatest shear in al] the oceans. These strong currents play a major role in 
oil spills and dispersion of other pollutants, and also transport significant amounts of 
sediment, having a dramatic effect on beach erosion. Additionally, the modeling of 
nearshore circulation assumes a special significance for military operations, where the 
success or failure of an amphibious assault, or a special forces operation is directly related 
to our ability to forecast the conditions within this unique operational theater. 

The physical understanding and mathematical modeling of hydrodynamic processes 
in the nearshore zone have substantially improved during the last two decades. Yet, most 
existing models are formulated based on a monochromatic wave description and have been 
tested and calibrated with laboratory data, due to difficulties associated with measurements 
within this harsh environment. In this study, monochromatic formulations are extended to 
random wave fields by using a probabilistic approach (Thornton and Guza, 1983). The 
numerical models developed here are compared with data acquired during the DUCK94 
(October 1994) experiment conducted on a barred ocean beach at the U.S. Army Corps of 
Engineers Field Research Facility near Duck, North Carolina. The objective of this 
dissertation i$ to improve our current understanding of the vertical and cross-shore 
distributions of both longshore and cross-shore currents, providing the first steps towards a 
comprehensive three-dimensional model. 


In Chapter I, observed vertical structure of the mean longshore current over a barred 


beach is examined to test the hypothesis that the turbulent bottom boundary layer of the mean 
longshore currents is logarithmic and to investigate the influence of wave-breaking generated 
turbulence on the mean longshore current profile. The relationship between bottom shear 
stress and bottom roughness, including the influence of ripples and mega-ripples is also 
investigated. This Chapter consists of a journal article accepted by the Journal of 
Geophysical Research. 

The vertical structure of the mean cross-shore current (undertow) over a barred beach 
is investigated in Chapter I. A surface mass flux model, including contributions from white 
water surface rollers on breaking wave crests (Lippmann and Thornton, 1997), is used to 
predict the cross-shore distribution of depth-averaged undertow. The impact of using linear 
and non-linear wave theories to estimate surface mass flux 1s investigated. The vertical 
structure of the undertow is modeled using an eddy viscosity closure to estimate the turbulent 
shear stress. The effects of different depth-dependent formulations for the eddy viscosity as 
well different boundary conditions on the vertical structure of the undertow are evaluated. 

In Chapter IV, the models of the vertical profiles of longshore and cross-shore mean 
currents developed in the preceding Chapters are combined to formulate a quasi three- 
dimensional model to describe the cross-shore distribution of the longshore current. This 
model includes turbulent mixing due to the cross-shore advection of mean momentum of the 
longshore current by the mean cross-shore current (Putrevu and Svendsen, 1994), and 
contributions from breaking wave rollers. The effect of introducing the observed cross-shore 
variation of bottom shear stress (Chapter I) on the longshore current prediction 1s 


investigated. 


Each of the three dissertation Chapters represents a separate paper for publication. 
As such, they are self-contained with their own introduction, theory description, data results, 


conclusions, and references. 
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II. VERTICAL PROFILES OF LONGSHORE CURRENTS AND RELATED BED 


SHEAR STRESS AND BOTTOM ROUGHNESS 


(This chapter consists of a journal article accepted by 


the Journal of Geophysical Research in May, 1997) 


VERTICAL PROFILES OF LONGSHORE CURRENTS AND RELATED BED 
SHEAR STRESS AND BOTTOM ROUGHNESS 
A.F. Garcez Faria, E.B. Thornton, T.P. Stanton, C.V. Soares’ and T.C. Lippmann’ 
Naval Postgraduate School, Monterey, California 93940 
ABSTRACT 

The vertical structure of the mean wave-driven longshore current over a barred beach 
is examined on three strong current days during the DUCK94 experiment, and it is found that 
the bottom boundary layer is well described by a logarithmic profile (mean correlation 
coefficient for all 22 profiles, 0.98). The logarithmic profile fits better in the trough where 
turbulent bottom boundary layer processes predominate, than over the bar where breaking- 
wave induced turbulence generated at the surface modifies the profile. The surface layer in 
the presence of waves is well described by adjusting the logarithmic profile for the 
intermittent presence of water and adding the alongshore component of the mass transport 
velocity (slope of the least-squares linear regression between model predictions and 
observations, 1.005 and rms error of 7 percent). 

Bed shear stresses calculated from logarithmic velocity profiles are equated to a 
quadratic bottom shear stress formulation. The associated bed shear stress coefficients vary 
by more than an order of magnitude across the surf zone (0.0006-0.012). Bottom roughness 
was measured throughout the nearshore using a sonic altimeter mounted on a moving 


platform. The bed shear stress coefficients are positively correlated with bottom roughness 
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(linear correlation coefficient, 0.6). A higher linear correlation coefficient (0.8) is obtained 
by subtracting skin friction from the total bed shear stress. 
INTRODUCTION 

Knowledge of the bottom boundary and surface layers is fundamental to 
understanding nearshore hydrodynamics and sediment processes. For steady flow, such as 
in a river, the bottom boundary layer is well described by a logarithmic profile. Only limited 
wave-driven, longshore current vertical profiles have been measured. Visser (1986) measured 
wave-driven longshore currents in a laboratory experiment using micro-propeller and laser 
Doppler velocimeters and found profiles approached a logarithmic form. In a similar 
laboratory wave-driven longshore current experiment using laser velocimeters, Simons et 
al.(1992) verified that the vertical profiles tended to be logarithmic. 

The velocity profile of a steady current is modified by the presence of waves. The 
Superposition of waves on the mean current produces enhanced bottom friction (e.g., Grant 
and Madsen, 1979; Christoffersen and Jonsson, 1985; Myrhaug and Slaattelid, 1989; and 
Sleath, 1990 ). As a result, the vertical gradient of the mean current near the bed is increased 
and a more uniform profile can be expected throughout most of the water column. The 
influence of waves is inversely proportional to water depth and wave frequency, with 
decreasing importance of nonlinear interactions between waves and currents within the 
wave-bottom boundary layer with increasing depth and frequency. 

Turbulence induced by breaking waves modifies the vertical profile of longshore 
currents. The downward momentum mixing produced from wave breaking-injected 


turbulence results in a more uniform velocity profile. Fredsoe and Deigaard (1992), and 


Church and Thornton (1993) argue this would result in an increased bed shear stress. 

In the following, mean longshore current profiles obtained over a barred beach are 
examined with the objectives of: |) testing the hypothesis that the turbulent bottom boundary 
layer of the mean longshore currents is logarithmic, 2) investigating the influence of the 
surface layer on the mean longshore current profile, and 3) examining the relationship 
between bottom shear stress and bottom roughness, including the influence of ripples and 
mega-ripples. 

THEORY 

The vertical profile of longshore currents is significantly affected by the bottom 
boundary and surface layers. The bottom boundary layer assumes more importance since it 
determines the general logarithmic profile shape over most of the water column. However, 
processes in the surface layer can modify the profile in the presence of waves and wind. 
Bottom Boundary Layer 

Neglecting molecular viscous stresses, the alongshore momentum equation ( y- 


direction) is written 


Opvy | Opvu | Opv? , Cpwy _ _op 


Ot ox Oy Oz Oy (1) 


The velocities are expanded into mean, turbulent, and wave-induced components, 
ae ee u. and w= w+ w, where (i = 1,2) refers to horizontal coordinates (x,y) and 
the mean vertical velocity 1s assumed equal to zero. After time averaging, (1) can be 


simplified with the aid of the following assumptions: 1) straight and parallel contours, e.g., 


=O) =0 (overbar indicates time averaging); 2) steady state conditions, 1.e., <C) ie 
y t 

3) wave-induced and turbulent velocity components are statistically independent 
(uncorrelated), i.e., 4 =0; 4) horizontal turbulent momentum flux is small compared to 
Opvu , Opvu 


wave induced momentum flux (Stive and Wind, 1982), i.e., , and can be 


x Ox 
neglected. Contributions due to momentum mixing caused by interactions between cross- 
opVU ites 
shore and longshore currents (a usually not negligible (Svendsen and Putrevu, 
x 
1994; and Garcez Faria et al., 1996). Nevertheless, they are neglected here for simplicity 
sake. Applying these assumptions, the alongshore momentum equations can be written 


_opwv _ dpvi | dpwy 


2 
OZ Ox Oz (2) 
which says that the sum of cross-shore changes in the wave-induced alongshore momentum 
flux and vertical changes in the wave-induced Reynolds stress are balanced by vertical 
changes in alongshore turbulent shear stress. 
Applying a first order turbulence closure, the alongshore turbulent shear stress can 


be defined 
mae OV 
42) = ~pw = pu (3) 
where y, is the turbulent eddy viscosity, which is assumed to be uniform with depth. 
Substituting (3) into (2) 


Ot,(Z) _ dpva | dpwr 


4 
Oz Ox Oz “) 





The cross-shore gradient of vi is constant in the alongshore direction for straight and 
parallel contours. Within the surf zone, the shallow water approximation holds, thus ti is 
assumed independent of depth. The wave-induced Reynolds stress term (wv) can arise from 
sloping bottom effects as well as from wave amplitude gradient effects and it would have a 
non-zero contribution even for linear wave theory (Deigaard and Fredsoe, 1989). 
Nevertheless, Rivero and Arcilla (1995), showed that wy is a linear function of depth, and 
hence its vertical gradient is a constant. Consequently, the nght hand side of (4) is 
independent of depth and the shear stress profile is determined by integrating (4) over depth, 


to give 
14) = [2,0)-t,(-A)] = (5) 


Linearly varying shear stress occurs in flows driven by uniform hydrostatic pressure 
gradients, such as in steady open channel flows, which are well described by a logarithmic 
velocity profile. Therefore it is hypothesized that a steady, uniform, turbulent boundary layer 
flow over a rough surface in the alongshore direction can be sectea by a logarithmic 


profile: 


V(z) = = In eae (6) 


oO 


where z iS positive upwards from the surface, A is the mean water depth, « is the Von 
Karman constant (0.4), v, is the alongshore shear stress velocity and z, is the physical 


roughness height, determined by bottom topography and sediment grain size. When waves 
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are present, nonlinear interactions between waves and currents within the bottom boundary 
layer increase the bottom shear stress. Following Grant and Madsen (1979), this additional 
stress can be modeled by an apparent roughness height z_, that is analogous to, but larger 
than z,. 

The mean bottom shear stress ( t,(-h)) is related to the longshore shear stress 


velocity (v,) through 


a) aa (7) 


In addition, a bed shear stress coefficient, Cr can be calculated assuming a quadratic bed 


shear stress relationship 





t(-h) = pc, Gee vy ey (8) 


and combining with (7), gives 


2 
Vv 


Cy = (9) 


(u2 + v2)'2y 


which includes contributions from both steady (U,V) and non-steady ( “7, V,u,v) velocity 
components. 
Surface Layer 

The surface layer is governed by waves and wind. Wave effects are manifested in 
three ways: 1) an undulating boundary is imposed on a mean current, which is being 


measured in an Eulerian frame; 2) the mass transport velocity in the crest-trough region due 


I 


to obliquely incident waves contributes to the longshore flow; 3) modifications of the profile 
associated with wave-breaking generated turbulence and the eventual contribution of wave 
rollers to the mass transport. 

The shape of the mean current profile in this layer is determined mostly by the fact 
that the current meter is intermittently in and out of the water in the wave crest-trough region, 
e.g., for a linear wave the current meter is out of the water half the time at the mean sea level 
(MSL hereafter), and the time-averaged current is only 50 percent of the expected value from 
the logarithmic profile. To account for this, the surface elevation probability distnbution 
function (pdf hereafter) is applied to the expected mean current profile in the absence of 
waves. The percent of time the current meter is in the water is given by 1-P(1), where P(n) 
is the cumulative surface elevation pdf. In an Eulerian frame of reference, the modified mean 


current within the surface layer is given by 


V(z) = [1-P(n)] V(z) (10) 


where V’(z) represents the logarithmic profile in the absence of waves (6). For moderate 

wave conditions in deep water, the surface elevation pdf is well described by the Gaussian 

pdf. As will be seen, the measured pdf in the surf zone are slightly positively skewed from 
the Gaussian distribution. 

The mass transport velocity associated with obliquely incident waves can contribute 

to the longshore flow in the upper layer. Assuming irrotational flow, the mass flux (transport) 

in the direction of wave transport, M, can be evaluated by considering separately two regions 


in an Eulerian frame of reference (Philips, 1977). In the first region, from the bottom to the 


[2 


MSL, the contribution is zero for irrotational flow. Within the second region, from the MSL 
to the water surface (1), # is not defined for linear wave theory. A Taylor series expansion 


about z=0 is used to extend defined values of “(0) at the surface, giving a second order 
approximation to M. The mass transport for a single wave is interpreted to this order, in an 
Eulerian reference frame, as due to a uniform velocity confined to the crest-trough region. 


Applying linear wave theory, a mass transport velocity is defined by 


M 
"pay a) 


where A is the wave amplitude. 

For random waves, the wave amplitudes can be described by the Rayleigh distribution 
to a first approximation, even within the surf zone (Thornton and Guza, 1983). The only 
waves that contribute at any elevation z, will have an amplitude A > Izl. Assuming 
directionally narrow banded waves, the ensemble averaged mass transport velocity profile 
in the direction of wave travel is obtained by applying the wave amplitude probability density 


function 


co 


<U(z)> = [ U(A)p(A)dA (12) 


Izl=A 


where < > represents ensemble averaging. The Rayleigh probability density function is given 


by 





- € ons where 0<A<h (13) 


and H, is the rms wave height. Substituting (13) back into (12), applying linear wave 


theory, and performing the integration, gives 


2!z! \2 
H (== 
O Ams .2Iz1 # VT 2Izl 
20>) = = —— eg ™ + 1. erfc(— or |Izli<h 14 
«) 8 tanh kh A 2 Tee ) J Oe 
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where erfc(x) is the complementary error function, @ is the wave radial frequency and k is 
the radial wavenumber. The alongshore component of the mass transport velocity is defined 


by 


<V(z)> = <U(z)>sinO (15) 


where 8 is the mean incident wave angle with respect to the beach normal, for the assumed 
directionally narrow banded waves. 

The mean longshore current within the surface layer, as observed in an Eulerian 
reference frame, 1s modeled by adding the corrected logarithmic profile for the cumulative 


surface elevation pdf (10) with the alongshore component of the mass transport velocity (15) 


Vidz) = L1-P(n)] Vz) + <U(e)>sinO (16) 


DUCK 94 EXPERIMENT 
The measurements described here are part of the comprehensive nearshore DUCK94 


experiment conducted during October 1994 at the U.S. Army Corps of Engineers Field 
Research Facility (FRF), Duck, North Carolina. The FRF is located on the Outer Banks, a 


barrier island formation with no major coastal structures to obstruct nearshore flows. The 
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beach is a two-bar system with a dynamic inner bar (30-120 m offshore) and a secondary bar 
with lower amplitude (300-400 m offshore). The mean foreshore slope of the beach is 
approximately 0.08 (1:12) and the slope offshore of the bars is approximately 0.006 (1:170) 
(Lippmann et al., 1993). The mean tidal range is 1.0 m. Sediments within the surf zone are 
well sorted with a mean grain size of 0.2 mm. Sediments on the foreshore are poorly sorted 
with larger mean grain size (>0.4 mm). 

The data analysis presented is for the October phase of the DUCK94 experiment. 
The weather during October was climatologically characterized by three distinct phases: 
weak currents and winds from north (4-9 Oct.), relatively strong currents from north (0.6- 
1.0 m/s) caused by a storm with predominant winds and waves from north (10-17 Oct.), and 
variable currents and winds from north/south (17-21 Oct.). For the first phase (4-9 Oct.), the 
currents were weak and barely above sensor accuracy. During the last phase (17-21 Oct.), 
currents and winds were highly variable and a hole developed in the bar associated with a np 
current system (observed with dye), such that the basic assumptions of steady state and 
Straight and parallel contours are violated. Within the second phase (10-17 Oct.), 
observations were limited due to the Coastal Research Amphibious Buggy (CRAB) not being 
able to operate during the height of the storm on the 15th and being limited to the area inside 
the bar on the 13th, 14th, 16th, and 17th due to large waves. Therefore, the data selected for 
analysis are from 10-12 Oct. (Fig. 1), during the strong longshore currents period when 
observations spanned the entire surf zone and conditions approximate the assumptions of 
steady state and straight and parallel contours. Within this period a logarithmic profile is well 


defined. 
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A specially designed sled was used as a platform to mount instruments (Fig. 2). The 
sled is constructed of a 3 x 4 m six-inch aluminum-pipe frame with two 5 m length, 20 cm 
diameter pipe runners. This low-profile structure was stabilized by 180 Kg of lead weight 
plus approximately 450 Kg of sand inside the runners. In addition, there were four fins (45 
cm wide) extending 60 cm into the sand to insure that the sled did not move while on station. 

Currents were measured using a vertical stack of 8 Marsh-McBimey two-component 
electromagnetic current meters (ems, hereafter) with 2.5 cm diameter spherical probes 
mounted along a 2.5 m mast (Fig. 2). The ems elevations above the bed were 23, 42, 68, 101, 
147, 179, 224 and 257 cm. The ems were displaced at least 1 m from the sled, and the sled 
was oriented such that the vertical stack of ems was on the “up-current” side of the sled to 
avoid flow contamination by the sled structure during observations. The ems were pre- and 
post-calibrated in a tow tank at the Naval Postgraduate School with an agreement of 1.9 
percent in gain. An in situ determination for the offset is used, which was obtained by 
reversing the orientation of the ems on a very slow longshore current day (8 Oct.) by turning 
the sled around and returning it to the same location (within 1 m) within one hour. The in situ 
determined offsets were within 1 cm/s. 

The sled orientation was determined using a digital compass mounted on the sled 
with accuracy O(1 degree). Measured two-component velocities were reduced to a shore 
normal right-handed coordinate system (positive offshore and to the south) by using compass 
data and adding at each sled position any deviation of the contour line, as measured by the 
CRAB, from a shore parallel direction. Velocity errors associated with the rotation of the 


coordinate system were smaller than the determined offsets, and therefore were neglected. 
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Waves and mean water level were measured using an array of five pressure sensors 
configured in a 3 m square with sensors at each corner and one in the center. 

The data were digitally encoded on the sled to 14 bit precision at 36 samples/second 
and transmitted to shore via a fibre-optic cable where signals were monitored and recorded. 
Short cables from the sensors to the data acquisition system on the sled (< 7 m) resulted in 
low noise ems and pressure sensor signals. An armored cable, married to the sled chain 
tether, provided power and controller signals for the instruments via two conductors and 
returned the digitized signals and video via a fibre-optic line. 

The sled was towed to the farthest offshore location for the first run (approximately 
160 m from the shoreline) by the 11 m high, motorized, three-wheel CRAB. A four-wheel 
drive forklift pulled the sled shoreward 10 to 30 meters for subsequent runs (each run was 
nominally one hour). Five to eight runs were made across a transect each day. The data were 
acquired during daylight to early night, which happened to span the high tide during this 
period. 

The morphology of the bottom (bathymetry) was measured at various scales from the 
CRAB. Large-scale variations of bathymetry were obtained by using an autotracking laser 
ranging system to measure the CRAB position approximately every meter with a vertical 
accuracy of less than 3 cm rms. Small-scale vertical bottom variations relative a the CRAB, 
including ripples and megaripples, were measured with a 1 MHZ sonic altimeter mounted 
on the CRAB, 70 cm from the bed. The altimeter has a 3.4 degree beam width which 
translates into an approximate 4 cm footprint and a nominal sampling rate of 25 Hz, which 


resulted in a sample spacing of 2-4 cm (dependent on CRAB speed) with mm vertical 
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resolution and accuracy less than 2 cm (Gallagher et al., 1996). The decrease in accuracy 
relative to resolution is due to the changing reflective surface owing to the bed dialating or 
sediment transported along the bed as waves pass overhead. The CRAB survey and altimeter 
measurements were combined to obtain a high resolution description of the bottom 
(Thornton et al., 1997). Contour plots of the bathymetry for the days selected for analysis 
show the alongshore contours in the vicinity of the sled measurements to be essentially 
Straight and parallel (Fig. 3; 11 Oct. is typical). Bathymetry for the three selected days 
(Figs. 4-6, upper panel) show a pronounced bar progressively moving offshore and 
significant small-scale morphology in the trough. Areal variations were determined using a 
500 KHz side-scan sonar also mounted on the CRAB. 

Meteorological information of wind, air temperature, atmospheric pressure and sea 
surface temperature were recorded simultaneously at the seaward end of the 600 m long FRF 
pier and atop the FRF building in front of the pier. 

DATA RESULTS 

The data are qualitatively sorted by location into the two regions of over the bar and 
in the trough. This sorting allows a better identification of the possible correlations among 
variables, as wave-breaking, which is a major controlling factor within the surf zone, 
significantly changes for these regions. 

Since the logarithmic velocity profile hypothesis is a bottom boundary layer concept, 
information from ems near the surface influenced by the effects of waves and wind as well 
as from coming in and out of the water are not included. A criterion is established such that 


only data from ems below (MSL- A) are considered to define the logarithmic profile of 
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the bottom boundary layer. This criterion assures that the ems used in the analysis came out 
of the water less than 0.25 percent of the time based on a Gaussian distribution, which is 
conservative for the measured positively skewed distributions. 


The rms wave height is approximated by H_ = ¥80° where ©” is the variance 


calculated from the surface elevation time series. Surface elevation was calculated by Fourier 
transforming a one-hour pressure record, applying a linear wave theory transfer function to 
the complex Fournier amplitudes in the frequency domain, and inverse transforming to obtain 
the surface elevation time series (Thornton and Guza, 1982). 

All 22 vertical profiles of longshore currents obtained during these three days are 
analyzed. The profiles are based on the measurements by four to seven ems over the vertical. 
The em closest to the sea bed was not used because of malfunction. Mean alongshore 
velocities are 1-hour averaged data with the exception of three runs, run 7 on 11 Oct. and 
runs 6 and 7 on 12 Oct., which are 40-minutes averaged data. In examining the bottom 
boundary and surface layers, the ems data are treated Separately. 

For the bottom boundary layer, logarithmic profiles are fit to the data based on a 
linear-regression least-squares method. The value of z, is calculated from the z intercept 
of the linear-regression on a semi-log plot of z versus V(z) , and the shear stress velocity v, 
is calculated from the slope. 

C f values are calculated using measured velocities (u,v) in (9) time-averaged over 
the record length, with v, determined by least-squares fit. C f values are dependent on the 
elevation of the measured u,v values. Typically, C f is calculated using measured u,v values 


referenced to one meter above the bed. A sensitivity analysis on the variation of C ¢ tO the 
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elevation of the selected gage was performed, showing that the variation of C, was 
dominated by the time-averaged alongshore velocity (V). This term was not only the largest 
in the denominator of (9), but also varied the most with depth. Therefore, to minimize the 


depth-dependence of C the depth-averaged velocity calculated from the logarithmic 





Zi 
velocity profile (V, = 7 - | J) is used to specify the dominant time-averaged 


alongshore velocity. The em isi at the elevation of 1 m above the bed 1s used to measure 
the smaller terms (total cross-shore velocity u and the non-steady component of the 
alongshore velocity )+v¥ = v-V), which have weak vertical variation outside the bottom 
boundary and surface layers. The total alongshore velocity is obtained by adding V, to the 
measured non-steady component V+¥v. The use of this method resulted in a mean variation 
of C, with depth of only 7 percent with a standard deviation of 6 percent and a maximum 
variation of 24 percent, provided the gage used to measure the cross-shore velocity and non- 
steady component of the alongshore velocity was not in the surface layer (z > MSL - H,,,). 
Calculated C (once and v, values are listed in Table 1. 


Error estimates of z, and v, based on the linear regression correlation coefficients 


(C) are calculated using (Gross and Nowell, 1983 and Cacchione et al., 1987) 


on 
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where 7 is the number of ems used for the regression, and t is the Student’s t 


(1-2, )-Gi2) 


distribution for (1-a@)confidence interval with (n-2)degrees of freedom. 
The uncertainties of C f values can be determined from the error estimates of v. and 


the time and depth-averaged velocity calculated from the measurements 
0 


ee al? 
Vm = - | [(v* +u7)\?y] dz (mean velocity, hereafter) 
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where v, is the relative error of the shear stress velocity at a given confidence level 


calculated from (17) and V__ is the relative error of the mean velocity determined by 


y= 1(cm/s) 


Morr 


, assuming that the mean velocity has a constant absolute error equal to the 
in situ oe offset (1 cm/s). Error estimates for Cr , Z,, and v, are listed in Table 2, 
for a 95 percent confidence interval. 

The measured cumulative surface elevation pdf is used to describe the surface layer, 
and compares reasonably well with the Gaussian cumulative pdf (Fig. 7). The Gaussian 
distribution has zero skewness and kurtosis equal to three. The measured skewness values 
ranged from 0.26-0.57 and kurtosis values ranged from 2.7-3.5, indicative of weakly 
nonlinear waves. 

In the surface layer, the upper ems sometimes came in and out of water, which 


causes noise in the ems outputs. To eliminate this noise, the current velocities were set to 


zero when the ems were within 5 cm of the surface as determined from the surface elevation 


Zh 


time series. 

The combined surface and bottom boundary layer solution is the logarithmic profile 
given by (6) from z, up to (MSL - H_,,), and the modified logarithmic profile by the 
observed cumulative surface elevation pdf plus the alongshore component of the mass 
transport velocity given by (16) above this level. As a typical example, the profile for the 
seventh run of 10 Oct. is enlarged in Fig. 8, where the solid line is the logarithmic profile and 
the dashed line the modified profile within the surface layer. 

The observed and model predicted velocity profiles at successive offshore positions 
(runs) that the sled occupied during a transect are shown in the upper panels of Figs. 4-6. The 
largest deviation occurs for the fourth run of 11 Oct., which generates an outlier in the 
parameter estimates 1f the data from the em at an elevation of 42 cm above the sea bed (Fig. 
5) is included in the linear regression. A comparative analysis for all ems for this run showed 
the energy density spectra to be nearly uniform (average variance of 0.072 m?and standard 
deviation of 0.013 m*) with exception of this em (variance of 0.021 m7). Thus, data from 
this em is disregarded and a new regression 1s calculated with the remaining six ems. This 
procedure eliminates the outlier. 

DISCUSSION 

Three days of the DUCK94 experiment are examined when strong longshore currents 
occurred. Mean longshore current profiles obtained using four to seven ems spaced from 42 
to 257 cm above the sea bottom are used to test the validity of the logarithmic profile 
hypothesis within the energetic surf zone region. 


Bottom Boundary Layer 


ae 


A high correlation coefficient for the linear regression, is commonly accepted as an 
indicator of the validity of the logarithmic approach (Grant et al., 1984; Gross et al., 1994; 
Li, 1994). The linear correlation coefficients for all profiles ranged from 0.95 to 0.99 
(Table 2), with an average value of 0.98, and the largest deviations occurring over the bar, 
where wave breaking was strong. Anomalous high correlation coefficients can be obtained 
for profiles calculated using a small number of ems for the linear regression, such as 
occurred for stations very close to the shoreline due to shallow water. Thus, a high 
correlation coefficient, although necessary, is not sufficient to validate the logarithmic 
profile approach. 

Other measures of how well the logarithmic model describes the mean alongshore 
current profile within the bottom boundary layer are the uncertainties of the calculated shear 
stress velocity v, (17) and apparent roughness height z, (18). These uncertainties reflect 
both the correlation coefficients and the number of ems used in the regression. An example 
is comparing the fourth run of 11 Oct. with the seventh run of 10 Oct. (Table 2) . The latter 
has a higher correlation coefficient (0.991), but only 4 ems were used in the linear regression; 
thus uncertainties in the calculated v, and z, values are larger for this run, even though the 
former has a lower correlation coefficient (0.987). Therefore, the uncertainties are not biased 
by an anomalous high correlation een due to a small number of ems used in the 
regression. The degenerate case is when only two ems are used for the linear regression 
(eighth run of 11 Oct.), resulting in a correlation coefficient of one, which does not allow 
the calculation of the uncertainties for v, and z,. Consequently, data from this run are 


disregarded. 
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There is some uncertainty of the exact distance of the ems from the bed, because the 
sled runners sank into the sand an unknown amount depending on the bearing capacity of the 
bed (theory and divers suggest 3-8 cm) and owing to the undulating bottom, particularly 
over mega-ripples. To test the sensitivity of the logarithmic profile due to the uncertainty in 
elevation, the elevation of the ems was shifted +5 cm in steps of 1 cm and the linear 
regression least-squares-fit recalculated until the best correlation between the data and the 
logarithmic model predictions was obtained for each run (Grant et al., 1984). 

As the bathymetry is not uniform over the transect covered by the sled, no unique 
zero-shift value should be expected to optimize the correlation coefficient for all runs. 
Another effect on the error calculation to be considered with an ad hoc decrease (increase) 
in the elevation of the ems is that an em previously located above (below) the model cut-off 
level (MSL - H,,) can be relocated to a new elevation below (above) this level. For 
example, an em previously neglected by the model is now included in the linear regression. 
The zero-shifting modifies the error estimates of v, , z,, and C, (17, 18, and 19) by 
changing the correlation coefficients and the number of ems used for the regression. 

With the uncertainties involved, the zero-shifting method was only applied whenever 
a significant reduction (greater than 10 percent) was obtained in the shear stress velocity 
error estimate. Applying this criteria, only three runs were shifted (by the same - 5 cm): 
seventh runs of 10 and 11 Oct., and third run of 12 Oct. The sensitivity of the ems to 
elevation from the bed as inferred from the zero-shifting method for these three runs is 
shown in Table 3. The general trend being that both the shear stress velocity and apparent 


roughness height errors decreased with increased negative shift, although the improvement 
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is only significant for the shear stress velocity. The values listed in Tables 1 and 2 for these 
runs were calculated after applying the zero-shift. 

The observed and model predicted velocity profiles at successive offshore positions 
(runs) that the sled occupied during a transect are shown in the upper panels of Figs. 4-6. The 
data agree well with the model indicating that the contribution from momentum mixing 
arising from interactions between cross-shore and longshore currents ( seVu ) neglected for 


simplicity sake, does not significantly modify the vertical profile of longshore currents, 





suggesting that is independent of depth. 


opVU 
Ox 

The largest discrepancies between measured and modeled profiles occur over the bar, 
where wave breaking is most intense (x-distance between 220 and 240 m). The surface layer 
during wave breaking is typified by an intense production of turbulence that eventually is 
dissipated in the shear layer at the lower boundary of the surface roller for spilling type 
breakers and more intense injection of turbulence for plunging type breakers. The increase 
of turbulent mixing due to wave breaking produces a more uniform vertical profile of the 
mean longshore current within the bottom boundary layer for a given bottom shear stress, 
compared with profiles in the absence of breaking. Therefore, larger discrepancies between 
observations and logarithmic profile predictions (lower correlation coefficients) would be 
expected for increased turbulent mixing caused by wave breaking. 

To test this hypothesis, percentages of waves breaking were determined from video 
recordings. The number of waves breaking is determined using the methods of Lippmann and 


Holman (1991), while the total number of waves is found applying the zero-up-crossing 


method to the surface elevation time series (See Lippmann and Thornton, 1997, for details). 
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The general decrease of the correlation coefficient between observations and logarithmic 
profile predictions with increasing percentage of wave breaking (Fig. 9) supports this 
conclusion. 
Surface Layer 

The mean current profile within the surface layer is modified in the presence of winds 
and waves. The winds act directly to generate longshore currents via the alongshore surface 
wind stress component, and indirectly by generating obliquely incident waves which then 
force the longshore currents via changes in the radiation stress. Whitford and Thornton 
(1993) measured the various terms in the alongshore momentum balance including wind 
stress and wave forces during the SUPERDUCK experiment (15-18 October 1986) at the 
same location as the DUCK94 expernnment. They found an average wind force to wave force 
ratio for this four-day period of 0.11 (range 0.02 - 0.33) for a mean alongshore component 
of wind speed of 3.6 m/s (range 1.9 - 5.3 m/s). During the three-day period studied here the 
mean alongshore wind speed was 9.1 m/s (range 7.7 - 12.6 m/s) resulting in an average wind 
force to wave force ratio of 0.21 (0.10 - 0.50). Despite the wave force dominance 
characterizing wave-driven longshore currents, the wind force contribution is not negligible. 

Winds and currents are approximately from the same direction (Fig. 1), thus the wind 
force effect 1s to increase the alongshore currents. As the logarithmic profiles are fit to the 
data based on a linear-regression least-squares method, the direct wind force effect is already 
included in the logarithmic profile. The indirect wind force effect of wave generation 
resulting in increased wave heights (higher H_,,. values) is also included in the alongshore 


component of mass transport velocity. The theoretical exponential decay of wind-induced 
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current speed with depth was not observed within the surface layer due to the dominant effect 
of the ems intermittently being out of the water, and therefore is not included in the model. 

The modification of the mean current profile within the surface layer is modeled here 
by correcting the logarithmic profile predictions for measurements in an Eulerian frame with 
an undulating boundary and adding the alongshore component of the mass transport 
velocity (16). As few ems were located above the (MSL - H_,,.) level for each run, a local 
correlation coefficient as an indicator of the validity of this approach is meaningless. 
Therefore, the validity of the model is evaluated by comparing model predictions with 
observations for the ensemble of ems located within the surface layer. 

The linearity of (16) allows an evaluation of the contribution of each term 
(intermittent-wetting and alongshore component of mass transport velocity) separately. For 
the intermittent-wetting term only, a plot of predicted versus observed velocities (Fig. 10) 
shows good agreement, but with the modified logarithmic profile by the cumulative surface 
elevation pdf slightly under-predicting the observed velocities. The rms error and slope of 
the least-squares linear regression between predictions and observations are respectively 13 
percent and 0.96, giving a slope-error of -4 percent. Next, including the alongshore 
component of the mass transport velocity (Fig. 11) results in reduced rms and slope errors 
of 7 percent and 0.5 percent. Despite the small magnitude of the mass transport term (on 
average only 11 percent of the intermittent-wetting term), its addition corrects the small 
under-prediction of using only (10) and improves the overall agreement with the data by 
reducing both the rms and slope errors. 


Although the results obtained by this first attempt (to the knowledge of the authors) 
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to model mean longshore currents within the surface layer are encouraging, improvements 
can be achieved by including contributions from at least two physical processes. The first, 
arises from the additional mass transport associated with the presence of wave-breaking 
generated rollers, and will increase the model-predicted velocity within the surface layer. The 
second is associated with the directional spreading of random waves and it will have an 
opposite effect, reducing the model-predicted velocity. Preliminary calculations indicate that 
these terms are of the same order of the mass transport term, and thus an order of magnitude 
smaller than the intermittent-wetting. The good agreement with the data obtained by applying 
this simple model is attributed to the dominance of the intermittent-wetting term and 
compensating effects of neglected contributions from the presence of wave rollers and 
directional spreading. 
Bed Shear Stress Coefficient 

The bed shear stress coefficient (C -) from the quadratic friction model is an 
important parameter in both nearshore hydrodynamics and sediment transport. The simplest 
formulation for longshore currents assumes steady state wave conditions and straight and 
parallel bottom contours, and results in an alongshore balance between cross-shore changes 
in wave-induced momentum (radiation stress, S,,) with the bottom shear stress, to give 


(Thornton, 1970) 
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where |ié| is the magnitude of the total velocity vector. The sediment transport formulation 
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of Bailard and Inman (1981), based on the work by Bagnold (1966), states that the immersed- 


weight sediment transport rate 1s proportional to the local rate of energy dissipation 


” . eC, |.) (K,, i K,) (21) 


where ee and e are dimensionless time-varying vectors associated with the bedload and 
suspended load transport rates. It is noted that both the longshore currents and sediment 
transport rates are direct functions of C,. Nevertheless, the physical processes governing 
these phenomena are inherently different. The longshore currents are mainly controlled by 
form drag of the bed forms and non-linear interactions between waves and mean current. 
Sediment transport is mainly related to skin friction due to sediment grains (Smith, 1977; 
Dyer, 1980), nIGuen the presence of waves increases the sediment entrainment rate (Grant 
and Madsen, 1979), thus increasing the transport rate. 

The bottom shear stress coefficient (C f ) varied by an order of magnitude across the 
surf zone, with the values offshore and over the bar in the order of 10°, while the values in 
the trough were in the order of 10°. An attempt was made to find empirical relationships 
between C, and measured physical parameters commonly used throughout the literature 
such as _ (ratio of near-bottom wave velocity magnitude from linear wave theory to 
mean current speed), the rms bottom roughness (r) and percent of wave breaking. 
Surprisingly, no statistically significant correlation was found between C, and sy 

C 7 Was found to be negatively correlated with percentage of wave breaking, with a 


linear correlation coefficient of -0.71, which is statistically significant at the 99 percent 


confidence level (Fig. 12). This is contrary to the theoretical analyses of Fredsoe and 
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Deigaard (1992) and Church and Thomton (1993). An effect associated with wave breaking 
in the surf zone is the generation of surface rollers that increases the mass transport within 
the surface layer, thus increasing the cross-shore return flow, or undertow, and consequently 
the average total velocity (u7+ v*)'“v. This effect reduces the C , Values calculated using (9) 
for increasing wave breaking. In addition, decreased C, with increased wave breaking may 
be physically related to the planing-off of wave-ripples due to increased near-bottom 
velocities associated with wave breaking and consequent reduction of form drag, as well as 
to stratification caused by an intense vertical gradient of suspended sediment near the sea 
bed. This stratification can diminish the turbulent flow intensity in the bottom boundary 
layer, leading to reduced bottom stress (Smith and McLean, 1977a; and Taylor and Dyer, 
1977). 

Cy would be expected to be related to bottom roughness due to enhanced form drag. 
The bottom roughness 1s examined by calculating wavenumber spectra of the bed. To 
calculate wavenumber spectra, the unevenly spaced data from the combined CRAB surveys 
and altimeter measurements are linearly interpolated to evenly spaced 2 cm increments of the 
cross-shore distance. The small-scale morphology in general shows large cross-shore 
variation; as a consequence, the condition of spatial homogeneity (stationarity) required for 
calculating averaged spectra is not met. Therefore, continuous bottom roughness 
wavenumber spectra are calculated for 20 m cross-shore segments at increments of 1 m 
across the surf zone. 

Lowest wave numbers are filtered by subtracting a third-order polynomial best-fit 


curve from each 20 m section. A 10 percent cosine-taper data window is applied to decrease 
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spectral leakage. The spectra are summed over three wavelength bands (0.2 - 1.0 m, 1.0 - 
3.33 m, and 3.33 - 20 m) plus the total band (0.2 - 20 m), resulting in 160, 28, 10, and 198 
degrees of freedom for each band respectively. The wavelength bands chosen are based on 
examination of individual spectra. The spectra were generally broad, indicating that several 
ripple wavelengths coexisted as a result of newly formed ripples combined possibly with 
residual ripples from the past to form a complex series of ripple patterns, plus the effects of 
alignment of ripples relative to the cross-shore measurement axis. The rms height of each 
band is calculated as the square root of the variance within each band. Note that variances 
sum, not rms heights, such that the rms height of the sum of the three bands is calculated 
from the square root of the sum of their variances. The rms bottom roughness (Figs. 4-6, 
lower panel) is dominated by longer wavelength features primarily associated with mega- 
nipples, with a small contribution from the superimposed shorter wavelength ripples O(1 m). 
The general trend being that bottom roughness was smoothest offshore and over the bar 
where wave ripples were planed-off due to higher near-bottom velocities, with increased 
roughness within the trough associated with mega-ripples (see Thornton, et al., 1997, for 
details). The rms bottom roughness (r) is calculated from the total band (0.2 - 20 m), 
and thus includes contributions from both large and small-scale morphology. Calculated 
values are listed in Table 1. | 
Roughness measurements were made in the cross-shore direction only. To relate C f° 
calculated from the alongshore bed shear stress, to roughness, it is assumed that the bed 
forms are quasi-isotropic, which is not unreasonable for mega-ripples in the trough. 


However, this assumption is violated for long-crested wave ripples, such as occurred for the 
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seventh run of 11 Oct. when large wave ripples (observed in the side-scan sonar) were 
oriented parallel to the flow (alongshore direction) which results in over-estimating the rms 
roughness associated with longshore currents at this cross-shore position. Altimeter data for 
12 Oct. is not available from 210-260 m due to the high false return levels induced by strong 
scattering of the acoustic pulses from air-bubbles entrained by waves breaking over the bar. 
Therefore, the CRAB survey data is used as the profile in this section, and hence the bottom 
is not as highly resolved in this area, compromising roughness estimates for the third, fourth, 
and fifth runs of 12 Oct. Thus, roughness data from these four mins, as well data from the 
eighth run of 11 Oct. (only two ems below the MSL - H_,, level) are disregarded. 

If data from runs with logarithmic profile correlation coefficients less than 0.98 with 
resulting relative errors for C 7 greater than 50 percent (Table 2) are also disregarded, only 
9 profiles could be used, making the ensemble too small to infer any statistically reliable 
relationship between Cr and bottom roughness. Therefore, it was decided to use data from 
the remaining 17 profiles with the associated 95 percent confidence error estimates to 
examine the relationship between C, and r. The bed shear stress coefficients are positively 
correlated with bottom roughness normalized by mean water depth (Fig. 13), with a linear 
correlation coefficient of 0.63, which is statistically significant at the 99 percent confidence 
level. 

The scatter of data observed in Fig. 13 may be due to correlating the non-synoptic 
velocity measurements with bottom roughness measured once in the morning, prior to the 
positioning of the sled for the first station. Wave forcing quantified by deep water wave 


height (H/_), period of peak wave frequency (7), and wave direction (0) changed during the 
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period of observation (Fig. |) due to both wind and tidal variations. Changes in wave forcing 
have a direct effect in the measured velocities and an indirect effect in bottom roughness due 
to modification of wave ripples associated with variations in wave height. 

For steady currents (open channel flows) C f depends only on bottom roughness when 
the bed is hydraulically rough, and several empirical relationships (power-laws) are available 
throughout the literature. One of the most widely used is the Manning-Strickler equation 


(Sleath, 1984, equation 5.86) 
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where f,,,, 1s the Darcy Weisbach friction coefficient and k, is the bed roughness length scale 
or equivalent Nikuradse roughness of the bed. 

In our definition of C f (9), contributions from both steady and non-steady velocity 
components were taken into account, whereas (22) applies to steady flow only. Nevertheless, 
during the period being analyzed the observed strong longshore current is not only the 
dominant term, but is also well represented by a logarithmic profile (6). The flow can be 
classified as hydraulically rough as calculated Reynolds numbers AX in which p = 
molecular kinematic viscosity coefficient) were greater than 70 for all runs (Schlichting, 
1979). Thus, if measured rms bottom roughness is representative of the equivalent Nikuradse 
roughness of the bed (k.=r), one might expect C, to be related to (=) by a power law 
similar to (22). A plot of predicted C, by the Manning-Strickler equation (22) 1s included in 


Fig. 13, which shows an order of magnitude agreement between the data and predictions. 
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These results suggest that the calculated rms bottom roughness (7) is representative of the 
bed roughness length scale or equivalent Nikuradse roughness of the bed (X_). 

The structure of the flow inside the wave boundary layer (WBL hereafter) is out of 
the scope of this paper as no measurements were made within this region. Outside the WBL, 
z, is the proper length scale and not z, (Grant and Madsen, 1986). Therefore, the equivalent 
Nikuradse roughness of the bed (k.) is replaced by an apparent bed roughness length scale 
(k_). The flow is hydraulically rough during the period being analyzed, and thus k, can be 
assumed to be given by ss = z, (Nikuradse, 1933). Assuming that k_ is representative of 
the bed roughness length scale when waves and currents are present, an empirical relation 
(power law) between C f and is obtained by applying a linear regression least-squares fit 


to the data to give 


aaa 
C, = 0.011 on aks (23) 


which is plotted in Fig. 14. The rms error between measured and predicted C f by (23) is 18 
percent, which is within the 95 percent confidence band of the measurements. This is not 
surprising as the coefficients of (23) were determined by fitting the data in a least-squares 
sense. Nevertheless, the correlation between measured C, and is 0.80 which is 
Statistically significant at the 99.5 percent confidence level indicating that a power law 
relation between these parameters is not fortuitous. 

An independent check of (23) is obtained by including the results from velocity 


profiles measured by Grant et al. (1984) on the northern California continental shelf in 
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90 m depth during the Coastal Ocean Dynamics Experiment (CODE) using four vertically 
stacked acoustic-travel-time current meters (BASS) mounted on tripods. A plot of 
measured C, versus = for the CODE data (listed in table 2 of their paper) is included in Fig. 
14. The agreement between C f predicted by (23) and measured for the CODE data set is even 
better (rms error of 3.5 percent) than for the DUCK 94 data set which was used to determine 
the coefficients of (23). 

These results indicate that a single roughness length scale can be used to characterize 
combined flows over a movable bed, which is a basic assumption for all existing theoretical 
models. This assumption was shown to be valid for waves, currents, and combined flows 
over a fixed rippled bed by Mathisen and Madsen (1996). 

Another reason for the scatter of data observed in Fig. 13 arises from C f values being 
calculated from the total bottom shear stress, which has stress contributions from skin 
friction related to sediment grains, waves-current nonlinear interactions within the bottom 
boundary layer as well as form drag related to bed forms. Smith and McLean (1977b) 
linearly partitioned the total bed shear stress into skin friction and form drag, and found good 
agreement with data from the Columbia River. Nelson and Smith (1989), Wiberg and Nelson 
(1992), and Li (1994) subsequently applied the linear partition concept to several flume 
experiments with good results. Extending this concept to the surf zone environment requires 
including an additional component due to non-linear interactions between waves and 
currents within the bottom boundary layer to the total bed shear stress. Assuming that the 
linear stress partition is valid within the surf zone, the skin friction contribution can be 


: a 
removed from the total bottom stress and a new bed shear stress coefficient Cy is defined 
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C= = (24) 
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wherev, and v_ are the form drag and waves-current interactions shear stress velocities. 
d we 


The relationship between G, and the total bed shear stress coefficient (C;) can be 


determined from (9) and (24) 


ve 
CG = U-() 1G, (25) 


where is is the skin friction related shear stress velocity. 

As skin friction was not measured during the DUCK94 experiment, an attempt is 
made to isolate its contribution from the total bottom shear stress by applying a stress 
partitioning model. The probabilistic approach used to quantify bottom roughness does not 
allow the adjustment of the two empirical coefficients C,, and a, necessary to apply the 
Smith and McLean (1977b) model. Therefore, the empirical relationships obtained by Li 
(1994) by applying linear regression to the mean flow laboratory data of Kapdasli and Dyer 
(1986), Paola (1983) , as well as his own data obtained from mobile sand ripples, are used 


to estimate skin friction shear stress velocity from the total shear stress velocity obtained 


from the logarithmic profile 
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where R is ripple height, that 1s assumed here to be equal to the measured rms bottom 
roughness (R =r). For -- > 6.86, the predicted Vs by (27) is greater than the total 
measured shear stress velocity, and it is assumed that Ue is equal to 99 percent of the total 
shear stress whenever this occurred. Calculated values of v, are listed in Table 1. 

The recalculated bed shear stress coefficients G, by (25) show a higher linear 
correlation coefficient (0.76) with bottom roughness normalized by mean water depth (Fig. 
15). Theoretically this is expected as the removal of the skin friction component from the 
total stress should enhance the form drag contribution, and consequently increase the 
correlation between bed shear stress and bottom roughness. The improved correlation 
obtained by applying Li (1994) empirical relationships indicates that the linear stress 
partition concept introduced by Smith and McLean (1977b) can be extended to the surf zone 
environment, although a more detailed data set is necessary to validate these expressions. 
Apparent Roughness Height 

The change of mean current velocity profile due to the superposition of waves has 
been theoretically investigated over the last two decades and several models have been 
presented (Lundgren, 1972; Smith, 1977; Bakker and Van Doom, 1978; Grant and Madsen, 


1979; Fredsoe, 1984; Christoffersen and Jonsson, 1985; Myrhaug and Slaattelid, 1989; 


Sleath, 1991). A general approach common to all these works 1s splitting the mean velocity 


Si 


profile in an inner region close to the bed affected by turbulence produced in the wave 
boundary layer, and an outer region above it, where the current is not directly affected by 
wave-current non-linear interactions so that the logarithmic profile (6) may be assumed to 
apply. The increased bed shear stress associated with wave-current interactions is modeled 
by replacing z, with z_in (6). 


a 


. a . 
The apparent roughness increase — is theoretically expected to depend on the 
= z 
: |i, | i A, 
relative current Strength ze the relative roughness — (ratio of the near-bottom wave 
5 
semi-orbital excursion from linear wave theory to equivalent Nikuradse roughness), and the 


. e . ° e ie 
angle between wave and current direction (). Several empirical relations for — based on 
2 


oO 


these parameters have been obtained from laboratory experiments. Nevertheless, the 
applicability of these empirical relations to field data has not been verified, as previous field 
experiments, although containing reliable data on z_, lack corresponding measurements of 
bottom roughness (Grant et al., 1984; Cacchione et al., 1987; Lambrakos et al., 1988; 
Slaattelid et al., 1990). 

Three widely applied empirical relationships for <a (Coffey and Nielsen, 1986; 
Sleath, 1991; and Van Rijn, 1993) are tested using DUCK94 mes as measurements of both z, 
and bottom roughness were made during this experiment. For hydraulically rough flows, the 
physical roughness height can be estimated from z, = = (Nikuradse, 1933), a once again 
the measured rms bottom roughness is used to represent the equivalent Nikuradse roughness 
of the bed (k_ = r). 

No statistically significant correlation was found between observed <a and the 


Z 
jt,| A ; 
parameters “ee = and @. Consequently, none of these empirical relations are able to 
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accurately predict the apparent roughness increase for this data. These results, together with 
the scatter of available laboratory data on <a (Nielsen, 1992, Fig. 1.5.13) and difficulties 
associated with bottom roughness and viscid measurements within the WBL in the field, 
indicate that much remains to be done before a reliable empirical or theoretical relation for 
the apparent roughness height increase experienced by the mean current profile in the 
presence of waves is obtained. 

CONCLUSIONS 

The vertical structure of mean longshore currents on a barred beach is well described 
by a logarithmic profile for the three strong longshore current days examined. This 
hypothesis works better in the trough where turbulent bottom boundary layer processes are 
more dominant than over the bar, where breaking-wave induced turbulence generated at the 
surface modifies the profile. 

The modification of the mean longshore current profile within the surface layer in the 
presence of winds and waves is modeled by correcting the logarithmic profile predictions for 
measurements in an Eulerian frame with an undulating boundary for all ems located above 
the (MSL - H.,,,) level. The addition of the alongshore component of the mass transport 
velocity corrects the under-prediction of the modified logarithmic profile and improves the 
overall agreement with the data. Therefore, it can be concluded that this simple model 
provides a first order approximation that is sufficiently accurate to predict mean longshore 
currents within the surface layer. 

The data indicates that wave breaking inside the surf zone decreases the bottom shear 


stress coefficient which is contrary to the theoretical analysis of Fredsoe and Deigaard 
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(1992), and Church and Thornton (1993). This may be physically related to wave-breaking 
generated surface rollers that increases the cross-shore return flow (undertow) and 
consequently the average total velocity, decreasing C,values calculated using (9); as well as 
to the planing-off of wave-ripples due to increased near-bottom velocities and consequent 
reduction of form drag, and to stratification caused by an intense vertical gradient of 
suspended sediment near the sea bed. 

The bed shear stress coefficient (C f ) varied by an order of magnitude across the surf 
zone (0.0006-0.012), with the values offshore and over the bar O(10°), while the values in 
the trough were O(107). Thus, longshore current and sediment transport models that assume 
G 7 10 be constant or mildly changing should be revised. 

C, was found to be directly proportional to bottom roughness, and hence, bottom 
roughness is an important parameter to characterize the bottom boundary layer. The 
empirical relationships obtained by Li (1994) were used to remove the skin friction 
contribution from the total bottom shear stress. The improved correlation between bottom 
shear stress and bottom roughness obtained, although not conclusive to validate these simple 
expressions, indicates that the linear stress partition concept introduced by Smith and 
McLean (1977b) can be extended to the surf zone environment. 


An empirical relation between C, and apparent roughness length scale of the bed 


k 
normalized by water depth a was obtained by applying a linear regression least-squares 
fit to the DUCK94 data. Good agreement was also found between this relation and CODE 
data (Grant et al., 1984), indicating that a single roughness length scale can be used to 


characterize combined wave-current flows over a movable bed. 
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Surprisingly, no statistically significant correlation was found between observed 


| zZ ju,| A 
apparent roughness increase (—) and the parameters aa —! and om. Consequently, none 


of the three most used empirical relations fora (Coffey and Nielsen, 1986; Sleath, 1991; 
and Van Rijn, 1993) were able to accurately _— the apparent roughness increase. These 
results indicate that despite improvements in our knowledge of turbulent boundary layers 
during the last two decades, we are still not able to estimate bottom roughness for field 
applications, which prevents the use of an empirical relation such as (23) or any other 
theoretical expression to predict friction factor for combined waves and current flows. 
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Figure 1. Climatology for the three days considered (10, 11 and 12 Oct.). Currents were 
measured in the middle of the trough. H is the significant wave height, T is the period of 


peak frequency, and Angle is the mean wave angle relative to beach normal (0). 
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Figure 2. Photograph of the Sled being pulled off the beach by the CRAB (upper right) 


during DUCK94 experiment. Vertical array of ems are mounted along Sled mast (center 


left). 
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Figure 5. Same as in Figure 4 but for 11 Oct. The anomalous measurement point indicated 


by + in the fourth run (upper panel) is not included in the linear regression. 
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Figure 6. Same as in Figure 4 but for 12 Oct. 
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Figure 7. Comparison between measured and Gaussian cumulative surface elevation pdf for 
all the runs during the three days considered (10, 11 and 12 Oct.). Mean skewness and 


kurtosis for the three days are indicated. 
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Figure 8. Vertical profile combining bottom boundary and surface layer models for the 


seventh run of 10 Oct. The upper three ems sometimes were out of the water. 
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Figure 9. Correlation coefficient between observations and logarithmic profile predictions 
versus percentage of wave-breaking calculated from video data. The line represents a linear 


regression with a regression coefficient of -0.94. 
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Figure 10. Predicted (equation 10) versus observed longshore velocities within the surface 


layer. The solid line represents perfect agreement (slope=1) and the dashed line represents 


a linear regression with a slope of 0.96. 
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Figure 11. Same as in Figure 10 but including the alongshore component of mass transport 


velocity (equation 16). The slope of the linear regression is 1.005. 
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Figure 12. Bed shear stress coefficient ( C f ) versus percentage of wave-breaking calculated 


from video data. The line represents a linear regression with a regression coefficient of -0.71. 
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Figure 13. Bed shear stress coefficient (C, ) versus rms bottom roughness normalized by 
mean water depth (>). The dashed line is the Manning-Strickler equation (22) and the solid 
line is obtained by adjusting the coefficients of (22) by applying a linear-regression least- 


squares method to the data. 
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Figure 14. Bed shear stress coefficient ( C ; ) Versus apparent bed roughness length scale 


normalized by mean water depth ce, The line represents predicted C, by equation 23. 
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Figure 15. Bed shear stress coefficient without skin friction contribution (er versus rms 
; Fr : . 
bottom roughness normalized by mean water depth — The line represents a linear 


regression with a regression coefficient of 0.76. 
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Table 1. Logarithmic Profile Fitting Results. 











Cross- rms rms Apparent Shear Bed Skin 
D R_ shore Depth Wave i Bottom Toughness = stress_~—s shear friction 
a u_ Position height y = roughness height velocity _ stress 
H shear 
y on h ne rms iE a va coeff. 

(m) = (m)— (m) (cm) (o21)) AG ean, aka 

velocity 

va 

(cm/s) 
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10 4 201 3.45 0.89 0.26 oe 0.6 7al 0.0048 3.9 
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Il 2 270 3.16 E27 0.40 1.8 02 4.7 0.0025 29 
IEE: 3 244 3.16 Le 0.35 1.3 0.9 2.6 0.0043 4.0 
ll 4 230 375 0.98 0.26 7.4 1.5 6.3 0.0063 3.0 
i 5 187 3.24 0.89 0.27 Oy) Ne 5.2 0.0053 2.4 
ll 6 el 545 0.84 0.28 * 9.0 oie 4.0 0.0065 Les 
Ine 7 157 2.50 0.76 0.30 16.0 0.5 25 0.0032 1.1 
ll 8 146 EGS 0.66 0.40 jg 3.2 oe 0.0048 1.1 
127 1 298 351 123 0.35 0.5 0.02 1.6 0.0009 1.0 
Zan 2 253 2.93 1.30 0.44 len 0.02 23 0.0011 1e5 
3 oe 2.98 1.21 0.41 ie 0.2 8.6 0.0022 Dee 
12 4 225 3.54 1.04 0.29 4.2 6.9 8.1 0.0121 5.0 
a 5 210 3.51 0.97 0.28 33 2.3 6.1 0.0069 a7 
VAs 188 2.99 0.90 0.30 6.3 1.8 4.7 0.0055 22 
i 7 172 2.74 0.84 O31 11.0 6.8 32 0.0073 iS 
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Table 2. 95 percent confidence bands on v_, z,, and C, . Run 8 of 11 Oct is not 


included as only two ems are used for profile fitting. 


Correlation 95% Confidence Band 
Number Student’s t coefficient 
of ems distribution between Shear Apparent Bed shear 
Day Run used for observations stress roughness stress 
the t (n-2,1-a/2) and log-profile velocity height Coefficient 
regression predictions Vv, a Cr 
n e 

(%) CSE) (%) 
10 l i 2,095 0.983 17 2 oo 
10 2 5 2.553 0.953 43 oe 99 
10 3 2 2553 0.952 44 7.4 102 
10 + 6 B32 0.992 14 1.9 27 
10 5 6 2.132 0.992 13 1.9 26 
10 6 5 2.559 0.999 5 1.3 6 
10 7 4 2.920 0.991 28 2 56 
11 ] 6 2.132 0.969 27 3.6 58 
11 2 5 2353 0.963 38 5.8 87 
11 3 5 2.393 0.975 31 4.1 68 
IL) 4 6 2132 0.987 17 2.4 35 
11 5 6 2eboz 0.991 15 2 28 
ll 6 5 2.353 0.998 9 1.5 14 
1] a 5 2553 0.996 12 1.8 21 
2 l 6 282 0.966 29 3.9 . “So 
12 2 + 2.920 0.968 53 10 128 
12 5 5 2353 0.989 20 25 41 
12 4 6 2lge 0.998 6 es 10 
Ihe: 5 6 oae2 0.997 8 Ness 14 
2 6 5 2353 O1995 14 1.9 25 
12 7 5 2.353 0.992 17 22 30 
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Table 3. Sensitivity of the ems to elevation from the bed as inferred from the zero-shifting 





method. 
Correlation 95% Confidence Band 
Number coefficient 
of ems between Shear Apparent 
Day Run usedfor  zero-shift observations stress roughness 
the and log- velocity height 
regression (cm) profile Va a 
predictions 
n C (%) (x/+) 

10 7 3 0 0.985 72 Za 
10 ic 3 -] 0.986 66 2.0 
10 iy 3 -2 0.986 60 2.0 
10 7 3 -3 0.986 55 2.0 
10 a 3 -4 0.987 49 fee 
10 7 4 -5 0.991 5) 1.8 
10 7 3 +5 0.984 107 22 
11 7 4 0 03 87 3.0 
11 7 + -] 0.994 81 Zo 
11 7 4 -2 0.994 75 29 
11 7 4 -3 0.994 69 2.8 
11 7 4 -4 0.994 64 2.8 
11 7 5 -5 0.996 ao 2a 
11 c 4 +5 0.992 124 a2 
IZ 3 4 0 0.987 37 4.1 
12 3 4 -] 0.988 34 4.1 
IZ 3 4 -2 0.988 31 4.0 
12 3 5 -3 0.988 20 Be 
12 3 5 -4 0.988 19 oof 
12: 3 5 -5 0.989 17 3.6 
IZ 3 4 +5 0.985 56 4.4 
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I. MEAN CROSS-SHORE CURRENTS OVER A BARRED BEACH 


(This chapter consists of work to be submitted 


to the Journal of Geophysical Research) 
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CROSS-SHORE MEAN FLOW OVER A BARRED BEACH 
A.F. Garcez Faria, E.B. Thornton, T.C. Lippmann’, T.P. Stanton 
Naval Postgraduate School, Monterey, California 93940 
ABSTRACT 

The spatial distribution of the mean cross-shore flow (undertow) over a barred beach 
is examined with field data obtained on three energetic wave days during the DUCK94 
experiment. The vertical structure of the undertow is modeled using a turbulent eddy 
viscosity closure and includes the important effects of wave breaking (described using the 
roller concept) and convective acceleration of the current. Other than a more realistic 
description of observed turbulence variations, a depth dependent eddy viscosity does not 
improve the agreement between predicted and observed undertow profiles. The effects of 
using different boundary conditions is investigated by extending the monochromatic 
formulations of Stive and Wind (1986) and Svendsen et al. (1987) to random waves by 
ensemble averaging over the wave height distribution. The contribution of breaking wave 
rollers to the surface mass flux can be of the same order of the contribution associated with 
the organized wave motion and should not be neglected. The largest discrepancies between 
model predictions and observations occur over the sand bar, where the observed strong 
undertow jet is not predicted by the model. It is postulated that the downward momentum 
mixing induced by wave breaking produces a constricted return flow layer that enhances 
undertow velocities, a process not accounted for in the model. Better agreement is obtained 
for this region with an empirical relation that correlates the surface mass flux with the mean 


rate of wave-energy dissipation per unit area. 


‘Scripps Institution of Oceanography, La Jolla, California 92903 
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INTRODUCTION 

The influence of the mean cross-shore return flow, or undertow, on sediment 
transport and beach profile evolution has long been believed to be important (e.g., Bagnold, 
1940). A local vertical imbalance between the wave setup pressure gradient, that is uniform 
with depth, and the depth-varying wave radiation stress is conceptually responsible for 
driving the undertow (Dyhr-Nielsen and S¢rensen, 1970). In the last two decades several 
theoretical models for the vertical structure of the undertow for two-dimensional beaches 
have been developed (e.g., Dally, 1980; Svendsen, 1984; Dally and Dean, 1984; Stive and 
Wind, 1986; Svendsen and Buhr Hansen, 1988; Deigaard et al., 1991; Stive and DeVriend, 
1994; Haines and Sallenger, 1994). All models use an eddy viscosity closure scheme, and 
solve for the depth dependent undertow by integrating the cross-shore momentum equation 
twice over depth, which requires two boundary conditions to evaluate the integration 
constants. 

There is a general consensus throughout the literature of using local conservation of 
mass over the vertical as one boundary condition. Commonly, the second boundary condition 
is either the stress at the trough level (Stive and Wind, 1986) or the no-slip condition at the 
bottom combined with the steady streaming generated by the bottom boundary layer 
(Svendsen, 1984). Despite significant physical differences, both approaches reduce to the 
same form between the trough level and the top of the bottom boundary layer within the surf 
zone, as contributions from steady streaming and bed shear stress are outweighed by mean 
water slope and wave forcing gradients. 


The wave-induced onshore mass flux in the region between the wave crest and trough 
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is critical to predicting the magnitude of the undertow, which 1s predicted heuristically by 
adding the contribution from breaking wave rollers to the mass transport given by an inviscid 
wave theory plus (Svendsen, 1984). Dally and Brown (1995) found better agreement for 
laboratory generated regular waves breaking over a planar beach using stream function theory 
(Dean, 1974) compared with linear wave theory (LWT, hereafter), which tended to over- 
predict the depth-averaged undertow. On the other hand, Masselink and Black (1995), using 
field measurements from experiments at two planar beaches, found good agreement using 
shallow water LWT. 

Observations of onshore mass transport by waves in the crest-trough region are scarce 
owing to difficulties associated with making measurements in this region. As aresult, the 
performance of surface mass flux models is usually evaluated by assuming that mass is 
conserved over the vertical and comparing model predictions with measured return mass 
flux. This procedure dictates a priori determination of the height of the cross-shore return 
flow layer. Direct determination of the height of this layer requires measurements of the 
vertical distribution of mean cross-shore flow with a vertical resolution that so far has only 
been achieved 1n laboratory experiments with monochromatic waves using Laser Doppler 
Velocimeters (e.g., Nadaoka et al., 1989). 

Most existing undertow models show good agreement with laboratory data for 
monochromatic waves breaking over planar beaches when the depth-averaged mean return 
flow is adjusted to fit the data (instead of using predicted mass flux) and the magnitudes of 
the two largest dynamical forcing terms (wave setup and radiation stress gradients) are 


determined from data. 
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Validation of these models with field data has been limited by the lack of data. Smith 
et al. (1992), compared a model based on Putrevu and Svendsen (1991) solution for the 
vertical structure of the undertow to data from the 1990 DELILAH experiment on a barred 
beach near Duck, North Carolina, and found large discrepancies over the bar, where the 
model under-predicted observed velocities. This strong “undertow jet” over the bar was also 
present during an earlier field experiment (Sallenger et al., 1983) at the same site. 

Within the surf zone, wave-breaking generated turbulence dominates bottom 
boundary layer processes, with the turbulent shear stress maximum at the surface decreasing 
to a minimum at or near the bottom. King and Kirby (1994), using LDV data from a wave 
flume experiment, showed that the primary turbulence generating mechanism in the surf zone 
is due to wave-breaking at the wave-roller interface and that turbulence intensities decrease 
with distance from the surface. Cox and Kobayashi (1997), using laboratory LDV 
measurements of regular waves spilling on a rough, plane slope, observed that the shear 
stress distribution within the surf zone varies linearly with depth until the top of the bottom 
boundary layer, and that the eddy viscosity is small near the trough level, increases to a 
maximum about one-third of the depth below the trough level, and then decreases toward the 
bottom. Svendsen (1984), investigated the effect of introducing an exponentially varying 
eddy viscosity (u,) to the predicted vertical profile of the undertow. In his formulation, two 
free parameters are necessary to define the magnitude and decay rate with depth of p_, and 
despite the more realistic depth-variation of ,, only marginal improvements are obtained 
for the vertical profile of the undertow (Svendsen and Buhr Hansen, 1988). 


In this paper, field observations of vertical profiles of the mean cross-shore current 
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obtained over a naturally barred beach are used to test various models. Wave height 
transformation is based on an energy flux balance including the effects of rollers to describe 
wave breaking and a probabilistic description of wave heights. The mass flux is investigated 
using both linear and higher-order wave theories, and including the contribution from surface 
rollers, to estimate the cross-shore distribution of the undertow. The influence of depth- 
dependent formulations for the eddy viscosity and different choices of boundary conditions 
on the vertical structure of the undertow is examined. 
THEORY 

In the following, a right-handed coordinate system with origin at the shoreline (x 
positive offshore) and z positive upward from the sea surface is used. Solutions for the 
surface mass flux, setup and the vertical profile of the undertow are described assuming 
Stationary wave conditions, straight and parallel depth contours, and random waves that are 
narrow-banded in both frequency and direction. 
Wave Transformation 

Wave transformation is determined from a probabilistic breaking wave model that 
includes roller energy gradients in the energy flux balance (Lippmann et al., 1996, LBT 
hereafter). This model assumes that wave heights both inside and outside the surf zone can 
be reasonably described by the Rayleigh distribution (Thornton and Guza, 1983) and has 
been shown to give accurate results for random waves breaking over both planar and barred 
beaches. The model has two free parameters: o, the mean angle of the wave-roller interface, 
and y, ameasure of breaking wave intensity y = a where is the rms wave height 


and h is the local depth. Although the model is insensitive to the interfacial angle (kept 


(2 


constant at o = 10°for all runs) based on the results obtained by LBT at the same location 
of the DUCK94 experiment, it is necessary to adjust y to fit the data. 

Wave properties are ensemble averaged by integrating wave heights through the 
assumed Rayleigh distribution, p(H). Following Thornton and Guza (1983), the fraction of 
waves that are breaking is found by integrating through the breaking wave distribution 
p,(H) = WA) p(H), where the weighting function, W(H), is given by Whitford (1988). 
Surface Mass Flux 

The conservation of mass for straight and parallel contours with the boundary 


condition of no flow through the beach is given by 


n 
[ elU@) +a) +u(z)]dz = 0 (1) 


ral 


where overbar indicates time averaging, 7 is surface elevation, is water density, and the 
onshore horizontal velocity has been partitioned into mean, wave and turbulence 
contributions. In an Eulerian reference frame and assuming irrotational flow below the 
trough level, there is a net shoreward mean mass transport by waves limited to an upper 


region between the crest and trough that is given by 
1) 


q, = [ pe iz)dz (2) 
" 


where subscripts (c, t) refer to (crest, trough). There is an additional contribution to the mean 


mass transport (per surface area) above the wave trough (surface layer) that arises from the 
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presence of turbulent wave rollers (Svendsen, 1984), that is defined by 


Cc 
= 90 A— 


where, is the roller density, Lis the length along the wave covered by the roller cross- 
sectional area A, and c is the velocity that the roller mass is advected landward, assumed to 
be given by the wave phase speed. The roller contribution (3) to the surface mass flux is 
based on the calculated cross-sectional area of the roller. Earlier models (Svendsen, 1984; 
Deigaard et al., 1991) assume that this area is proportional to the rms wave height. As a 
consequence, the largest mean return flow would be predicted to occur at the breaker point, 
which is contrary to laboratory observations (Nadaoka and Kondoh, 1982). Dally and Brown 
(1995) solve for the area of the roller for the case of monochromatic waves breaking over a 
planar beach by simultaneously solving the depth-integrated and time-averaged energy, 
continuity, and cross-shore momentum equations and found good agreement with existing 
laboratory data. However, their model requires observations of the cross-shore distribution 
of either the mean return flow or setup to sancti the model, and therefore cannot be 
applied to the present data. 

Another recent roller model (Lippmann and Thornton, 1997; LT97 hereafter) presents 
an independent method for calculating qg_, and is calibrated with video observations of the 
cross-shore variation of the fraction of waves that are breaking. The LT97 model is based on 
the energy flux balance and describes energy dissipation following Deigaard (1993). Two 


free parameters in the model, B, the vertical fraction of wave height covered by the roller, 
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and , a measure of the average wave face angle, are adjusted to give a rms best fit to 
breaking observations. 
The onshore mass transport in the upper region is balanced locally by a mean return 


flow below the trough (undertow) 


N: 
dy +9, = - f PU@)dz = -pUgh, > 
-h 


where U, is the depth-averaged return flow and A, is the depth below the trough. 
Undertow 

The time-averaged cross-shore momentum equation, neglecting molecular viscous 
stresses, and for the assumed straight and parallel] contours and stationary wave conditions, 
is given by 


dpu? A Opuw _ _ op 
Ox OZ Ox 


(5) 
The horizontal and vertical velocities (u,w) are expanded into mean, turbulent, and wave- 
induced components, v= U+u+u and w= w+w, where the mean vertical velocity is 
assumed equal to zero. The time-averaged pressure, p=p[g(n-z)- Ww? - w?] , 1s obtained 
from the depth integrated vertical momentum equation, after neglecting contributions from 
the cross-shore gradient of the vertically integrated wave and turbulent shear stress (Stive and 
Wind, 1982). After substituting for the time-averaged pressure, (5) can be further simplified 


for the region between the top of the bottom boundary layer and trough level (middle layer), 


TS 


with the aid of the following assumptions: 1) wave and turbulent velocity components are 

statistically independent; 2) turbulence is near isotropic (Stive and Wind, 1982); 3) the wave 

RA eel 2 aaah iveroand rele loosyranalaaiieet 
Oz OK 


order eddy viscosity closure for the turbulent shear stress - puiw = p Ls . Applying these 
z 


shear Stress 1s given by 


assumptions, (5) reduces to 


= 772 
— [a*-w*] + pg— + ——— F(x) (6) 


where p,, is the time-invariant turbulent eddy viscosity. The forcing, F(x), in (6) is due to the 
cross-shore gradients of radiation stress, setup (setdown), nN, and convective acceleration of 
the depth-averaged undertow respectively. Here, F(x) is assumed to be independent of depth 
based on empirical evidence from laboratory studies (Nadaoka and Kondoh, 1982; Stive and 
Wind, 1982 and 1986). 

A solution for the vertical distribution of the mean undertow can be determined by 


integrating (6) twice over depth, to give 





d 


U(z) = F(x) i a er) i 

PH, p 
where C(x) and C,(x) are spatially varying integration constants. Both, the Stive and Wind 
(1986) and Svendsen et al. (1987) solutions are used, with modifications to their original 
developments to improve the derivation. In both these original works, a different solution for 


the undertow forcing (F) is used as they assume that the wave contribution to stress 1s small 
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compared with the Reynolds stress (paw « pw) resulting in an over estimation of the 
forcing by 52 (a2 - 2] (Rivero and Arcilla, 1995). Here, (6) will be used instead. 
Bs 
Although different formulations for the eddy viscosity variation with depth are 


investigated (Appendix), the simplest solution for a depth-independent eddy viscosity (u) 


is outlined here. For this case, (7) simplifies to: 


1 
pu 





= F 2 4 am, + 
U(z) 5 (x) z i Coe * Ce (8) 


First, following Stive and Wind (1986), C,(x) is determined by integrating (6) once 
over depth and solving for the shear stress at the trough level, and C,(x) by applying 


conservation of mass over the vertical (see Appendix for details) to give 


Ue) =U, + Pe) the + eye een ty 
z) = —————— x) i— + Za + — at = — 
“pu 2 6 pn 2 


where ce is the mean cross-shore bed shear stress that is calculated using a quadratic 


formulation ( oe = Cr u, U,, where C = 0.01 is aconstant friction factor and u, is the 


amplitude of the near-bottom wave-induced velocity). This monochromatic solution is 


quadratic in z with coefficients independent of wave height A, = = A, = h, and 
3h? - h? 


oF —— which can be extended to random waves by ensemble averaging over the 


wave distribution 
z 


~(z +h - A (10) 
pu 2 





<U(2)> = Up + (Ay2? +A,2 +g) — f Fp dH + 
pu 4 


Ts 


In the cross-shore direction, Te is small compared with wave-breaking dissipation (Thomton 
and Guza, 1983; Svendsen, 1984; Dally and Brown, 1995), and hence the vertical structure 
of the undertow is mostly determined by the coefficients A,, A,and A). The amount of 
curvature in the profile is a function of both F (calculated from wave quantities) and the eddy 
viscosity, .. Large values of F produce more vertical shear resulting in a parabolic profile, 
whereas large values of pp reduce the vertical shear producing a more uniform velocity 
profile with depth. 

The influence of the boundary condition choice to the vertical structure of the mean 
undertow is investigated by comparing the results obtained by (10) with the Svendsen et al. 
(1987) model that uses a no-slip condition at the bottom to replace the stress at the trough 
level as the second boundary condition. This no-slip condition is obtained by coupling the 
undertow model with a bottom boundary layer model (BBL, hereafter). Within the BBL, the 
flow is a combination of the steady streaming induced by the oscillatory motion and the 
undertow above the BBL, which results in a mean velocity at the top of this layer ( U, ) that 
is obtained by requiring continuity in velocity and shear stress between these two regions 


to give (see Svendsen et al., 1987 for details) 


, = eae 
<UG@)> = U, + (A,z?+A,z+A,)— ff FQ)p(HdH + BE (z +h) aD) 
PH PH 


oe a ee 
where A,= — is another coefficient independent of wave height and 1, is related to the 


ao . 
yeti: mon ~~ l aa a7 . 
steady streaming inthe BBL(t,. = - 5 Ou,, | ——., where f is wave frequency and 1, 
Ma Veh 


is the eddy viscosity inside the BBL). Their model further assumes that p,, is much smaller 
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than the eddy viscosity in the middle layer, which was recently verified in a laboratory 
experiment using LDV (Cox and Kobayashi, 1997). Following Putrevu and Svendsen (1993), 
the eddy viscosity inside the BBL is estimated by a UEZ6 A (ay 
hgh : 
Setup 
The setup gradient is a dominant driving force for the undertow. The setup is 


calculated by depth-integrating the time-averaged cross-shore momentum equation (5) from 


the bottom to the mean water level, to give 











0S. 0M Pe h) 0U, a h) on 7 
+ + e = - +h)—+ 
pvieee ea ae ae oo) 


where S$, is the wave radiation stress (Longuet-Higgins and Stewart, 1964), M_ = cq, is 
the momentum flux associated with wave rollers, and the last term on the lhs is the 
convective acceleration of the depth-averaged undertow. In the derivation of (12), it is 
assumed that the mean bed shear stress (Longuet-Higgins and Stewart, 1964) and the cross- 
shore wind stress (owing to small observed onshore winds in the data described later) are 
negligible. 
DATA 

Field measurements were acquired as part of the DUCK94 experiment (e.g. Garcez 

Faria et al., 1997 [Chapter 1]) conducted at the U.S. Army Corps of Engineers Field 


Research Facility (FRF) in Duck, North Carolina. The data selected for analysis are from 10- 


12 Oct. when strong cross-shore currents (0.05-0.4 m/s) caused by a storm with predominant 


winds and waves from the north were present. 


i 


The vertical structure of the current was measured with a vertical stack of seven two- 
component Marsh-McBimey electromagnetic current meters (ems, hereafter) mounted on a 
mobile sled at elevations 0.4, 0.7, 1.0, 1.5, 1.8, 2.2, and 2.6 m above the bed. The ems were 
horizontally displaced at least 1 m from the sled, and oriented such that the vertical stack was 
in the up-drift direction of the longshore current to minimize interference by the sled 
structure. The em offsets were determined in situ and found to be accurate within 1 cm/s. 

The sled orientation was determined using a digital compass mounted on the sled 
with accuracy O(1 degree). Measured two-component velocities were rotated to a shore 
normal right-handed coordinate system by using compass data and adding at each sled 
position any deviation of the bottom contour line from a shore parallel direction. A precise 
determination of the rotation angle is important to avoid contamination of cross-shore 
velocities by the observed strong longshore currents O(1 m/s). Velocity errors associated 
with the rotation of the coordinate system were comparable to the determined offsets, and 
therefore are neglected. 

For the first run on each day, the sled was towed by the Coastal Research Amphibious 
Buggy (CRAB) to its farthest offshore location seaward of the bar (approximately 160 m 
from the shoreline). A forklift on the beach pulled the sled shoreward 10 to 30 meters for 
subsequent measurement runs that are referred to in the text by sequential numbers within 
each day. Each data run was nominally one hour, and seven to eight runs were made across 
a transect during each day spanning the high tide during this period. 

Waves and mean water level were measured using an array of five pressure sensors 


mounted on the sled. Directional wave spectra were also acquired using a linear array of 10 
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pressure sensors in 8 m depth. Additionally, a 13 element cross-shore array of pressure 
sensors was used to measure wave heights spanning the width of the surf zone (Elgar et al., 
1997). The fixed array was located approximately 25 m to the north of the sled transect. 
These data were continuously sampled at 2 Hz. Video observations were used to measure the 
fraction of wave breaking along the same transect using the method of Lippmann and 
Holman (1991). 

Meteorological information of wind and atmospheric pressure were recorded 
simultaneously at the seaward end of the 600 m long FRF pier and atop the FRF building in 
front of the pier. The bathymetry was measured daily using the CRAB and it was found that 
the depth contours were nearly straight and parallel for the three days under consideration. 

RESULTS AND DISCUSSION 

Measured and modeled undertow flows are maximum on top and on the shoreward 
slope of the bar, with increasing magnitude from 10 to 12 Oct. as wave forcing increased 
with the approach of a storm (Fig. 1). The vertical structure over this region is the classic 
parabolic shape associated with strong wave-breaking turbulence (Svendsen, 1984; Stive and 
Wind, 1986). In the inner trough, the return flow 1s weak, and almost no vertical structure 
is seen. At the seaward slope of the bar, observed profiles on 10 Oct. (Fig. la) are nearly 
uniform with depth. It is noted there is significant bottom roughness in the trough of the 
barred profiles related to the development of lunate and long-crested mega-ripples as a result 
of the strong longshore currents (Thornton et al., 1997). The bar migrated off-shore during 
this period O(20 m), which was found by Gallagher et al. (1997) to be associated with the 


strong undertow. 


8] 


Model/data comparisons are evaluated by calculating absolute and relative root- 


mean-square (rms) errors. Absolute error has dimensional units and is defined by 


(13) 





where P. is model prediction, O, is the observed quantity, and N is the number of 


observations. The relative or percent error is calculated by 


E (14) 


rel 





and weights the difference between model predictions and observations against a measure 
of expectation that is represented by the observation. This statistic is not well-behaved for 
small O., which sometimes is the case for observed mean cross-shore velocities. Thus 
measurements smaller than the in situ determined offsets (+/- 1 cm s"') are excluded from 
(14). 

In the following, the accuracy of the wave transformation model is described first, 
which is then used to predict the surface mass flux, setup and vertical profiles of the 
undertow including contributions from surface rollers. 

Wave Transformation 

The rms wave height is approximated by H_ = 80° where o° is the variance 

calculated from the surface elevation time series. Surface elevation was calculated by Fourier 


transforming a one-hour pressure record, applying a LWT transfer function to the complex 
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Fourier amplitudes in the frequency domain, band-pass filtering by zeroing coefficients 
outside the range of interest (0.05Hz <f < 0.5 Hz), and inverse transforming to obtain the 
surface elevation time series (Thornton and Guza, 1982). 

The sensitivity of the models for surface mass flux, setup, and undertow to errors 
associated with the use of a random wave transformation model (LBT) is investigated by 
comparing with results using measured wave heights interpolated with cubic splines. The 
difference between mode] output and data interpolation methods is shown in Fig. 2 for the 
first run of 1] Oct., which has the largest rms relative error (7 percent). 

A plot of measured versus modeled H_ | for the three days analyzed (10-12 Oct.) 
shows good agreement (Fig. 3). The mean ¢_,, is 5 percent for all runs, with largest errors at 
any cross-shore position within 12 percent. Best fit values for y and €_, are summarized in 
Table 1. In general the transformation model well represents measured wave heights. 
Surface Mass Flux 
Estimating Surface Mass Flux from velocity measurements 

The surface mass flux is inferred from the measured return flow, which is sensitive 
to the depth of the return flow layer. For field experiments as well as for modeling efforts, 
the height of the return flow layer is assumed to be coincident with the trough level, h,. Here, 
h, is calculated using both linear and non-linear wave theories and compared to surface 
elevation estimates. Laitone’s second order approximation (Laitone, 1961) to cnoidal wave 
theory (Korteweg and de Vries, 1895) is used as this theory provides an analytical solution 
to h,, and the conditions for all runs being examined are within the range of validity of this 
H/(2*h) 

kh) 


theory determined by the Ursell number, > 0.6, where k is wavenumber (Laitone, 
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1962). The use of a non linear wave theory generally gives higher trough levels than LWT, 
and consequently smaller mean return flows ( U, ) for the same onshore mass transport (Fig. 
4). However, the effect on the calculated total return mass flux 1s minimal as an increase in 
h, is counterbalanced by a decrease in U,. No statistically significant difference at the 97 
percent confidence level is found between the results obtained using cnoidal and linear wave 
theory for al] 22 runs (Fig. 5). 

A problem associated with practical applications of cnoidal wave theory is related 
to its mathematical complexity, as it contains both the Jacobian elliptic functions and the 
complete elliptic integrals of the first and second kinds. An attractive alternative is the use 
of its asymptotic approximation, or hyperbolic wave theory (Iwagaki, 1968), which is valid 
when the complete integral of the second kind (K) is greater than 3, and the Ursell number 
is greater than 0.6 (Svendsen, 1974), which is the case for the data examined here. The 
calculated difference between observed surface mass flux given by LWT and Iwagaki’s 
approximation is €,,, = 7.9 percent, indicating that the hyperbolic wave theory is accurate 
enough to replace the more complex cnoidal wave theory. Applying Iwagaki’s solution, an 
expression for /, similar to LWT is derived that is a function of only local depth and rms 
wave height. 
eerie 
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The observed mean trough level is obtained by applying a zero-up-cross method to 


the measured surface elevation time series to define individual waves and then averaging the 
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minima between each up-cross interval. A comparison with measured mean trough level 
(Fig. 6), shows that this simple expression (15) provides a better estimate for the mean 
trough level within the surf zone (rms over-prediction of 4 cm) than LWT (rms under- 
prediction of 17 cm). It can also be noticed in Fig. 6 that discrepancies between predicted and 
observed mean trough level increase with rms wave height for LWT, while (14) provides an 
estimate of h, that does not appear sensitive to H__. 

Surface Mass Flux model 

The inferred surface mass flux obtained by numerically integrating the observed 
cross-shore mean currents from the bottom to the trough level given by (14) are compared 
to model estimates (4). g, 1s determined from the calibrated LT97 model and gq, is 
determined first using linear wave theory and then with non linear theory for comparison. 
Predicted g, are dependent on the percentage of breaking waves in the LT97 model. For the 
22 profiles examined, best fit values of B range 0.60-0.95 and wrange 1.1-1.9 (Table 1), and 
the error between observed and modeled percentage of breakers is €,,= 5.6 percent. An 
example of observed and predicted percentage of breakers is shown in Fig. 7 for the third run 
of 12 Oct., which has the largest rms error (9 percent). 

Inferred and predicted surface mass flux using linear wave theory are compared for 
the entire ensemble in Fig. 8. Due to significant changes associated with wave-breaking 
characteristics, the data have been divided into four regions: shoreward slope of the bar, 
seaward slope of the bar, trough, and foreshore. The rms relative error between observations 
and predictions using only q_, given by linear theory is 40 percent. Including the mass flux 


contribution from wave rollers, g_, improves the overall agreement (7ms error = 28 percent). 
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Largest values of q, occur at the shoreward slope of the bar (Fig. 8, lower panel), where 
wave-breaking is most intense. In this region, g, 1s on average 72 percent of gq, , but can be 
as large as 144 percent (fourth run of 12 Oct.). This is in accordance with earlier model 
studies (Svendsen, 1984; Dally and Brown, 1995), but contrary to Masselink and Black 
(1995) who contend that the roller contribution to the mass transport is of secondary 
importance based on the results from field experiments at two near planar beaches. 

The impact of using a non linear wave theory to calculate gis examined next. For 
this investigation, stream function theory (SFT, hereafter) 1s used as it provides a higher 
order solution for g, than cnoidal wave theory, and it has already been shown to give 
accurate results with laboratory data (Dally and Brown, 1995). Surface mass fluxes were 
calculated using SFT both including and excluding q,. For the field data examined, there is 
little difference between the use of LWT or SFT, with LWT giving values on average 8 
percent larger than SFT (Fig. 9). 

Errors in surface mass flux predictions could also arise from the misfit in estimating 
wave heights with a transformation model. Calculated rms relative error in mass flux 
between using LBT model predicted H_ _ and applying a cubic spline to the measurements 
for all runs is 9 percent, with a maximum relative error of 19.5 percent (third run of 12 Oct.). 
These results suggest that the surface mass flux is not overly sensitive to the choice of the 
random wave transformation model. 

Undertow Jets 
The onshore transport is observed to extend below the mean through level during 


high wave conditions and where a high percentage of wave-breaking occurs, such as over 
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of the bar (Fig. 4b), which violates a basic assumption of existing surface mass flux models. 
Inferred downward momentum mixing produced from breaking-injected turbulence is 
probably responsible for this anomalous onshore flow in the upper layer penetrating beneath 
the trough level (Masselink and Black, 1995; Garcez Faria et al. 1997 [Chaper 1]). To 
conserve mass, undertow velocities in this constricted return flow layer are enhanced 
resulting in an undertow jet. This phenomenon was qualitatively observed by Nadaoka et al. 
(1989) in a laboratory experiment using flow visualization techniques and a fibre-optic LDV 
system. In their experiment, obliquely descending eddies behind the breaking wave crest 
were observed to substantially transport vorticity into an otherwise nearly irrotational 
velocity field, resulting in an increased mass transport. 

On the shoreward slope of the bar where strong wave-breaking was observed (Fig. 
7), the model considerably under-predicts mass transport (Fig. 8). The effect of the 


de> | 4 acuegel 


constricted return flow layer is estimated for this region by . Enhanced 


Fobs 
undertow velocities account for an average of 39 percent of the iid Surface mass flux, 
and are in qualitative agreement with Nadaoka et al. (1989). 

These observed strong jets are not believed to be associated with three-dimensional 
circulation cells, as measured bathymetry was essentially uniform alongshore during the 
period being analyzed and no qualitative evidence of stationary np currents was observed in 
the video. 

Empirical estimate for Surface Mass Flux 
It was found empirically that the inferred surface mass flux (Q = pU,, h,) is 


significantly correlated at the 99.5 percent confidence level (linear correlation of 0.89 and 
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€_,,= 27 percent) with the mean rate of wave-energy dissipation per unit area, D (Fig. 10) 
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The spatially varying D is estimated from the energy balance equation using the LT97 model. 
This empirical relation provides a simple alternative to estimate the surface mass flux within 
the surf zone that easily can be incorporated into nearshore current models for field 
applications with the same accuracy of existing surface mass flux models. A better agreement 
is obtained for the region of the shoreward slope of the bar with this relation (Fig. 8-10). 
Setup 

Setup is calculated using a finite centered difference method to numerically solve 
(12), with the condition that the setup is assumed zero at the most offshore grid point. 
Contributions from each term of (12) are examined for the fourth run of 12 Oct. (Fig. 11), 
which corresponds to the most energetic period. The cross-shore gradient of the momentum 
flux associated with wave rollers is calculated using the LT97 model, the radiation stress 
term is calculated using linear wave theory (Longuet-Higgins and Stewart, 1964), and the 
convective acceleration of the current (2) is estimated by applying cubic splines to 
measured depth-integrated cross-shore mean return flow. The largest contributions are due 
to the radiation stress and roller terms. The convective acceleration term, although generally 
an order of magnitude smaller than the other terms, cannot be neglected, as it can be of the 


same order as the sum of the larger terms. 


The effects on setup by contributions from the roller momentum flux and using non 
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linear wave theory to calculate the radiation stress are examined with one example from 
each day (Fig. 12). The effect of the roller is to redistribute momentum laterally into the 
trough of the bar, shifting the point where the setup begins onshore. This shift has a 
significant impact on the setup/setdown profile within the surf zone. The use of a non linear 
wave theory (SFT is used here) to estimate the wave radiation stress does not appear to alter 
the setup profile dramatically (Fig. 12). 
Undertow 
Cross-shore variation of the eddy viscosity 

An implicit assumption of the closure model used in the undertow solution is that the 
eddy viscosity coefficient (1) is proportional to turbulence intensity. Conceptually, large 
cross-shore variation of this parameter is expected throughout the surf zone over a natural 
barred beach associated with significant changes in wave-breaking generated turbulence. 
Haines and Sallenger (1994), in an earlier field experiment limited to 3 ems over the vertical 
at this same beach, used a different model for the vertical variation of the undertow and 
found that best fit up for each location varied by more than an order of magnitude across the 
surf zone (u = 0.0055 - 0.075 m’ s"') in qualitative agreement with theoretical expectations. 
Smith et al. (1992) on the other hand applied a constant wu across the surf zone (u= 0.05 m’ 
s') and found good agreement with the DELILAH data. 

In order to investigate the cross-shore variation of the eddy viscosity, best fit up (Table 
2) were calculated using (10) by minimizing both the €,,. and €, at each cross-shore 


position occupied by the sled (Fig. 1). For all 22 runs, €,,. = 2.2cms” (€_,, = 19 percent) 


rel 


with maximum €,, = 6.1 cm s'(e_, = 44 percent) for the fourth run of 12 Oct. Observed 


rel 
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mean undertow velocities (U_, _) and splined rms wave heights were used in the calculations 


obs 
to isolate effects due to the choice of .. Using modeled rms wave heights (LBT96), €,,. = 
2.5 cms" (€_,, = 22 percent), which indicates that the vertical profile of the undertow is not 
sensitive to errors associated with the use of a random wave transformation model. Using the 
surface mass flux model to estimate U, results in rms errors of 5.1 cms” and 48 percent, 
which suggests that largest errors in the vertical profile of mean undertow are related to 
failure of existing models to correctly predict the surface mass flux in breaking waves. 
Models for the vertical structure of the undertow that use an eddy viscosity closure 
to solve for the turbulent shear stress (Reynolds stress), independent of the choice of 


boundary conditions and depth-dependence of the eddy viscosity, result in a general solution 


of the form: 


U(z) = U,, + a St) (17) 


where Use is a reference velocity that could be either the depth-averaged undertow or the 
velocity at the top of the bottom boundary layer. A nearly uniform vertical profile of 
undertow, as observed in the inner trough and seaward slope of the bar (Fig. 1), requires 
either no forcing (fF = Q) or an infinite eddy viscosity to be correctly reproduced. An 
example of the spatial distribution of each term and the total F(x) for the first run of 10 Oct. 
is shown in Fig. 13. The dominance of the setup gradient term is evident, and hence the 
model predicts offshore directed flow throughout the entire surf zone, which agrees with 
observations. It can also be noticed that modeled total forcing, although small offshore of the 


bar and within the trough where the undertow was near uniform, was never nil within the 
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cross-shore transect occupied by the sled, which was also the case for all the other runs 
during the period being analyzed. Minimum rms errors over these regions are obtained for 
large eddy viscosity values (Table 2), with a maximum cutoff level arbitrarily set at 
pp = 0.5. These unrealistically large values of wu are associated with the sensitivity of the 
eddy viscosity model to small errors in the dynamical forcing, which is calculated as the 
difference between two large numbers (setup and radiation stress gradients). 

The parameterization for the eddy viscosity is of the form u = 1 C, where / andC 
are characteristic length and velocity scales (Battjes, 1975). Several dimensionally consistent 
parameters for w are found in the literature, such as ys [Thornton, 1970]; 
h = [Battjes, 1975]; h Vgh [Stive and Wind, 1986]; an oy [Haines and 
Sallenger, 1984]. An attempt was made to relate best fit with these parameters after 
excluding from the ensemble runs that reached the cutoff level (Table 2). No statistically 
significant correlation at the 95 percent confidence level was obtained between these 
parameters and best fit u. On the bar crest and shoreward slope of the bar the mean undertow 
profile assumes a parabolic shape and best fit eddy viscosity approaches a constant value 
(u = 0.04). Although contrary to theoretical expectations, these results suggest that the use 
of a constant eddy viscosity across the surf zone provides a reasonable approximation for the 
modeling of the vertical structure of the undertow. 

Vertical variation of the eddy viscosity 

Here, the impact of both a depth varying p._ and different boundary conditions on the 
vertical structure of the undertow under field conditions are investigated. It is expected from 
measurements (King and Kirby, 1994; Cox and Kobayashi, 1997) that yp. should increase 
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from the bottom towards the surface with increasing levels of wave-breaking generated 
turbulence. Several mathematical formulations for the vertical variations of , were 
investigated and the three solutions that give best agreement with data are described in the 
Appendix. A comparison between the solutions for a linear (A6) and quadratic (A11) p, that 
give the smallest overall rms errors for the entire ensemble, together with the undertow 
solution (10) for a constant eddy viscosity with depth are shown in Fig. 14 for the two 
stations within each day that have largest observed vertical structure of the undertow. No 
significant improvement in the rms relative or absolute errors between observations and the 
solutions using variable , with depth is observed. These results confirm earlier modeling 
findings that a depth-dependent eddy viscosity does not substantially improve the description 
of the vertical structure of the mean undertow (Svendsen and Buhr Hansen, 1988; Nadaoka 
et al., 1989). 

The influence of the boundary condition choice to the vertical structure of the mean 
undertow in the middle layer is investigated next. Comparison between predictions by (10) 
that uses conservation of mass over the vertical and the stress at the trough level as boundary 
conditions and (11) that uses a no-slip condition at the bottom to replace the stress at the 
trough level as the second boundary condition are shown in Fig. 15 for the same stations used 
in Fig. 14. Again, no significant improvement in the total rms errors between observations 
and predictions by these two solutions 1s obtained, although there are noticeable differences 
in the predicted currents by each model as a function of depth. Equation (10) shows better 
overall comparison with data, while (11) better represents the structure of the flow in the 


lower half of the water column. Although models coupling the middle layer and BBL flows 
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(11) are expected to provide a more realistic description of the undertow structure close to 
the bed (Battjes et al., 1990), their applicability to field conditions is still limited by the lack 
of data to constrain their free parameters. 

CONCLUSIONS 

The predicted spatial distribution of mean cross-shore currents (undertow) over a 
barred beach are compared with field observations during the DUCK94 experiment to 
quantify the relative importance of contributions from the various terms in the cross-shore 
momentum equation and to identify physical mechanisms not yet incorporated in existing 
nearshore models. 

The surface mass flux model based on conservation of mass over the vertical and 
including contributions associated with wave rollers is compared with field observations 
obtained by integrating the measured return flow over depth. It was found that the roller 
contributions to the mass transport can be larger than the contributions from the organized 
wave motion when high energetic waves are present. An increase of 8 percent in the surface 
mass flux prediction is obtained using linear wave theory combined with contributions from 
wave rollers compared to the solution given by non linear streamfunction wave theory. 

The largest discrepancies between model predictions and observations of undertow 
velocities occur over the sand bar, associated with the failure of existing models to correctly 
predict the surface mass flux under breaking waves. It is postulated that the downward 
mixing produced from breaking-injected turbulence penetrating beneath the trough level 
results in a constricted return flow layer and enhanced undertow velocities (undertow jet), 


a process not accounted for in the model. Good agreement with data is obtained for this 
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region with an empirical relation that correlates the surface mass flux with the mean rate of 
wave-energy dissipation per unit area. This empirical relation provides a simple alternative 
to estimate the surface mass flux within the surf zone that can be easily incorporated into 
nearshore current models for field applications with the same accuracy of existing surface 
mass flux models. 

Setup (setdown) is a dominant forcing mechanism for the undertow and is calculated 
from the depth-integrated and time-averaged cross-shore momentum balance. It is found that 
despite the dominance of the setup (setdown) and radiation stress gradients in this balance, 
the contribution from the convective acceleration of the mean current is significant during 
energetic wave events. It is also shown that the inclusion of contributions from wave rollers 
results in an onshore shift of the point where the setup begins, which has a significant impact 
on the dynamical balance of forces within the surf zone. 

The vertical structure of the undertow is modeled using a turbulent eddy viscosity 
closure, and it is found that a depth-dependent eddy viscosity does not substantially improve 
the agreement with data under field conditions. The effects of using different boundary 
conditions to solve for the vertical structure of the mean undertow is investi gated by 
extending the monochromatic formulations of Stive and Wind (1986) and Svendsen et al. 
(1987) to random waves by ensemble averaging over the wave distribution (Equations 10, 
and 11, respectively). No significant improvement in the total rms errors between 
observations and predictions given by these two solutions is obtained. Despite the better 
representation of the undertow structure close to the bed obtained by using (11), which 


couples the middle layer and bottom boundary layer flows, its applicability to field 
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conditions 1s still limited by the lack of data to constrain their free parameters. 

In the inner trough and seaward slope of the bar, the measurements show almost no 
vertical structure for the mean undertow, which would require a local balance between setup 
gradient and wave forcing (F = 0) or a large eddy viscosity (1 - ©) to be properly modeled. 
Although small, model predicted forcing was never nil within these regions, so that 
unrealistically large values of » are required to mode] the observed uniform profile of the 
undertow. These unrealistically large values of pare associated with the sensitivity of the 
eddy viscosity model to small errors in the dynamical forcing, that is calculated as the 
difference between two large numbers (setup and radiation stress gradients). 

On the bar crest and shoreward slope of the bar, where strong wave-breaking was 
observed, the mean undertow profile assumes a parabolic shape and best fit eddy viscosity 
approaches a constant value (u = 0.04). These results suggest that the use of a constant 
eddy viscosity across the surf zone provides a reasonable approximation for the modeling of 
the vertical structure of the undertow. 

The absence of both setup measurements and observations of currents close to the bed 
are major limitations of the present study. Setup is an integral measure of the dynamic 
response of the mean water level to cross-shore gradients of momentum fluxes, and its 
gradient is an important driving force for the undertow within the surf zone. The lack of 
current measurements close to the bed prevents a quantitative evaluation of the effect of 
applying different boundary conditions to solve for the vertical profile of the undertow. 


APPENDIX 


In this appendix, the solution for the case of an eddy viscosity that varies linearly with 


v5 


depth, as suggested by Okayasu ef al. (1988) is derived in more detail. Following Stive and 
Wind (1986), the stress at the trough level and conservation of mass over the vertical are 
used as boundary conditions. The final solutions for other depth-dependent formulations for 
the eddy viscosities investigated in this study are also included. 


1. Linear: p, = « + Bz’ 


_ (Z+h) 
t 
primes are omitted hereafter for simplicity of notation. Assuming the undertow forcing to be 


In this solution, a normalized vertical coordinate is introduced z’ , and the 
constant over depth and using the eddy viscosity closure for the turbulent shear stress, a 
solution for the vertical distribution of the mean undertow is obtained by integrating (6) over 


depth, to give 


0 U(z) 
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where C, is an integration constant that is solved by evaluating (A1) at the bottom (z = 0) 


Cant. A2 


Integrating Al a second time and using A2, gives 


U(z) = Je) [z ~ G ince + Bz)] + ay inte sea a) are A3 


pp 


where C, is an integration constant that is solved by applying conservation of mass over the 


vertical, to give 
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substituting A4 in A3, the final solution for the vertical profile of the undertow is given by 
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2. Quadratic: pW, = a + B(z)'? 
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This parabolic form for the vertical variation of the eddy viscosity was suggested by 


Roelvink and Reniers (1994), and observed by Cox and Kobayashi (1997) with LDV in a 


laboratory experiment on regular waves spilling on a rough, plane slope. 
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In the above listed solutions, the eddy viscosity (y_) is a function of two non- 
dimensional parameters (wand B). The value of the eddy viscosity at the bottom ( Es 
given by a, which in this study is assumed to be represented by the molecular viscosity, due 
to the lack of reliable measurements of currents close to the bed under field conditions to 
better constraint this value. Although, the solution within the middle layer does not 
significantly change if the assumed constant eddy viscosity within the bottom boundary layer 
(Putrevu and Svendsen, 1993) is used for a. 
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Figure 1. Measured ( * ) and rms best fit model predicted (Equation 10) vertical profiles of 


mean cross-shore return currents (heavy line) superposed on bottom profile with meantrough 


level indicated by (0) for the entire period being examined. 
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Figure 2. LBT “worst case” model prediction of H,. versus cross-shore distance (dash- 
dotted line) compared with cubic spline interpolation (solid line) of observations ( * ), first 


run of 11 Oct. (upper panel). The bottom profile is shown in the lower panel. 
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Figure 3. Predicted (LBT model) versus observed H_.. The solid line represents perfect 


agreement and the dashed lines the +10 percent error bounds. 
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Figure 4. Observed vertical profiles of mean cross-shore currents ( * ) for the eighth run of 
11 Oct. (a) and for the fourth run of 12 Oct. (b). Horizontal heavy lines represent mean 
trough level (h,) and vertical heavy lines the depth averaged undertow ( U_) calculated using 


linear (dashed) and cnoidal (solid) wave theories. 
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Figure 5. Comparison between inferred surface mass flux (p U_, h,) from mean velocity 


measurements using linear and cnoidal wave theories to estimate the mean trough level, h,. 
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Figure 6. (Predicted - Observed) mean trough level, h,, versus observed H_ .. Positive 


values represent over prediction. 
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Figure 7. Lippmann and Thornton (1997) model predictions (solid line) and observed ( 0 ) 
percentage of waves breaking versus cross-shore distance for the third run of 12 Oct. (upper 


panel). The bottom profile 1s shown in the lower panel. 
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Figure 8. Inferred (p U 


0. 


,,2,) versus linear wave theory predicted surface mass flux (¢,,, 
only) for the entire ensemble of 22 runs (upper panel). The lower panel shows the effect of 


including contributions from wave rollers (q, + q,). 
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Figure 9. Comparison between predicted surface mass flux (¢, + q¢,) given by linear and 
streamfunction (Dean, 1974) wave theories. The dashed line represents perfect agreement, 


and the solid line represents a linear regression with a slope of 0.92. 


113 


SURFACE MASS FLUX [ Kg/ m/sec ] 
700 


seaward Slope of the Bar 


600 Shoreward Slope of the Bar 


Tr 7 
Trough 

F500 , 

= careeh 

— oresnore 

“7 400 

> 

3 

= 300 

O 

£ 

al 

a 200 

Lu 


100 





0 200 400 600 
INFERRED (SLED DATA) 


Figure 10. Inferred versus empirically estimated (Equation 16) surface mass flux. The line 


represents perfect agreement. 
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Figure 11. Terms of the momentum balance equation (12) versus cross-shore distance for the 


fourth run of 12 Oct. (upper panel). The bottom profile is shown in the lower panel. 
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Figure 12. Setup calculated using linear wave theory (dashed line), streamfunction theory 
with roller (dash-dotted line) and linear wave theory with roller (solid line) versus cross- 


shore distance for the fourth run of each day. 
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Figure 13. Undertow dynamical forcing terms (6) versus cross-shore distance for the first 
run of 10 Oct. (upper panel). The bottom profile is shown in the lower panel. The vertical 


dashed line indicates the cross-shore position occupied by the sled. 
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Figure 14. Comparison between measured undertow (*) and predictions using different 
formulations for the vertical variation of the eddy viscosity (constant = solid line; linear = 
dash-dotted line; and parabolic = dashed line) for the third and fourth runs of Oct. 10 (a, and 


b), Oct. 11 (c, and d), and Oct. 12 (e, and f). 
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Figure 15. Same as Figure 14, but comparing predictions given by Equation 10 (solid line) 


and Equation 11 (dashed line). 
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Table 1. Best fit model parameters 


Model: 


Run 
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0.32 
WZ 
HESZ 
0.32 
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0.34 
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Table 2. Wave conditions and undertow measurements 


depth- rms best fit 
rms wave peak mean tide averaged eddy 
Day Run height frequency direction undertow viscosity 
H, fp Qo Vobs LU 

(m) (Hz) () (m) (ms”) (m*s”) 

10 l 1.70 0.171 44 0.40 0.06 0.500 
10 Z 1.81 0.162 38 0.79 0.08 0.500 
10 3 1.81 0.162 38 0.98 0.16 0.058 
10 4 1.68 Oal5Z 36 0.77 0.14 0.017 
10 5 1.68 Or52 36 0.44 0.05 0.500 
10 6 1.44 0.162 24 0.04 0.05 0.500 
10 7 1.44 0.162 24 -0.12 0.05 0.500 
1] ] ZAlel 0.142 18 0.22 0.10 0.240 
I 2 1.88 0.142 16 0.64 0.12 0.250 
let 3 eS 0.142 1 0.86 0.19 0.033 
1] 4 1.70 0.142 18 97 0.19 0.049 
Ie 5 ee 0.142 18 0.85 0.05 0.063 
It 6 1.70 0.142 18 0.62 0.08 0.500 
1] 7 1.60 0.142 16 0.32 0.07 0.500 
1] 8 1.60 0.142 16 0.03 0.15 0.140 
|Z ] 1.9] 0.162 18 0.05 0.13 0.130 
12 2 2.29 0.142 2 0.34 0.20 0.240 
12 3 29 0.142 2: 0.71 0.30 0.032 
2 4 252 0.142 10 0.83 0.30 0.036 
12 » Zoe 0.142 10 0.84 0.20 0.035 
12 6 235 0.142 10 0.69 0.07 0.020 
|Z A eS 0.142 10 0.42 0.10 0.500 


H,, fp, and O&)were measures in 8m depth (FRF linear array) 
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A QUASI-3D MODEL FOR LONGSHORE CURRENTS 
A.F. Garcez Faria, E.B. Thornton, T.C. Lippmann’, T.P. Stanton, R.T.Guza’, and 
S. Elgar’ 
Naval Postgraduate School, Monterey, California 93940 
ABSTRACT 
The spatial distribution of mean longshore currents over a barred beach is examined 
by comparing predictions given by a quasi three-dimension model presented here with data 
acquired during the nearshore dynamics experiment DUCK94 at Duck, North Carolina. The 
model includes forcing due to breaking waves described using the roller concept (Lippmann 
and Thornton, 1997), alongshore wind stress, cross-shore advection of mean momentum of 
the alongshore current, and a full non linear bottom shear stress with a variable bed shear 
stress coefficient, C ir constrained by observations (Garcez Faria et al., 1997a [Chapter 1]}). 
Contributions from the alongshore wind stress are mostly evident offshore and over the inner 
trough of the sand bar due to the relative increase in the wind force to wave force ratio as 
wave forcing decreases over these regions. The advection of the momentum of the longshore 
current by mean cross-shore currents 1s shown to improve the agreement with observations 
within the surf zone, O[10 percent]. The use of a non linear bed shear stress formulation with 
a variable C f is shown to improve model/data comparison, O[20 percent], compared to the 
use of a constant C;. The largest overall improvement with observations is obtained by 


incorporating the roller contribution, O[50 percent]. 


‘Scripps Institution of Oceanography, La Jolla, California 92903 
“Washington State University, Pullman, WA 99164-2752 


124 


INTRODUCTION 

The longshore current maximum observed in the trough of a barred beach during the 
nearshore dynamics experiments DELILAH and DUCK94 at Duck, North Carolina, is not 
predicted by present theory. The simplest longshore current models balance cross-shore 
changes in the alongshore wave momentum (radiation stress) with the alongshore bottom 
shear stress. Waves break over the bar reform in the trough and again break on the foreshore. 
Wave breaking results in changes in the radiation stress predicting two jets, one over the bar 
and the other at the foreshore, which does not agree with the observed current maximum in 
the trough. An example of the measured and modeled wave height and longshore current 
distributions using DUCK94 data are shown in Fig. 1. The predictions suggest that a transfer 
of momentum 1s required to account for the current deficit in the trough. 

A number of mechanisms have been proposed to mix momentum laterally into the 
trough region to drive the longshore current. Traditional turbulent mixing, usually 
parameterized using classical eddy viscosity concepts associated with the shear of the 
longshore current (e.g., Bowen, 1969; Longuett-Higgins, 1970; Thornton, 1970), would 
require up gradient mixing of the longshore current which is not feasible. Battjes (1975) 
formulated turbulent mixing induced by breaking waves, where the scale of the turbulent 
mixing is the same order of the wave height or the local depth, which are much too short 
length scales to explain the observations. Smith et al (1993) described wave breaking as 
rollers that propagate with the wave at the phase speed; they applied a turbulent kinetic 
energy equation and argued that turbulence was diffused downward into the water column 


generating an additional alongshore thrust. They applied their formulation to the DELILAH 


ID 


data and partially explained the momentum deficit in the trough. 

Changes in the bottom shear stress due to turbulence being injected from the surface 
by breaking waves and modifying the vertical profile of the longshore current, which in turn 
modify the magnitude of the longshore current, were investigated by Church and Thornton 
(1993). They allow a spatially variable bed shear stress coefficient dependent on the breaking 
wave induced near bottom turbulence levels. ‘The model predicted cross-shore profiles of the 
longshore current improved agreement with observation compared with treatments using 
constant bed shear stress values, but did not completely account for the momentum deficit 
in the trough. 

Instabilities of the longshore current have been identified as a mechanism for 
turbulent mixing of the longshore current originally suggested by Bowen and Homan, (1989). 
Dodd and Thornton (1990) showed that if shear instabilities exist, there is an accompanying 
cross-Shore mixing of momentum. Putrevu and Svendsen (1992) carried out a numerical 
study of shear instabilities over various topography and using an order of magnitude analysis 
concluded that even a weak shear in the longshore current might be capable of producing 
significant mixing. Significant shear instabilities of the longshore current have been observed 
in the field (Oltman-Shay et al, 1989). Dodd et al. (1997) calibrated the amplitude of the 
shear instabilities using field measurements, and then calculated the Reynolds’ stress 
associated with the instabilities. They found that the mixing predicted due to shear 
instabilities to be in qualitative agreement with that required for modeled longshore current 
profiles to agree with observed profiles. 


Putrevu et al. (1995), using an order of magnitude analysis, suggested that alongshore 
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pressure gradients induced by alongshore variation in bottom topography could contribute 
substantially to the forcing for longshore currents. Reniers et al. (1997) examined this effect 
using DELILAH data and found that even a relatively small alongshore pressure gradient 
acting in the same direction as the wave forcing over the trough of the bar, can have a 
significant effect on the longshore current distribution. 

The advection of the momentum of the longshore current by mean cross-shore 
currents as a source of momentum mixing was suggested by Putrevu and Svendsen (1993). 
The longshore current is strongest toward the surface and decreases to zero at the bottom. 
The cross-shore mean current has an onshore transport in the wave crest/trough region and 
an offshore transport below (undertow). Svendsen and Putrevu (1994) showed that the net 
interaction could induce significant mixing of the momentum of the mean currents. 

The objective of this paper is to formulate a simplified quasi-three-dimension model 
to describe longshore currents assuming that the bottom profile has straight and parallel 
contours. Included in the formulation are forcing due to breaking waves described by a roller 
(Lippmann and Thornton, 1997) and alongshore wind shear stress, cross-shore advection of 
mean momentum of the alongshore current by the shear of the mean cross-shore current as 
suggested by Svendsen and Putrevu (1994) and a full non linear bottom shear stress. The 
observations acquired during DUCK94 will be used to test the model and to quantify the 
relative importance of contributions from the various terms of the alongshore momentum 
equation within the surf zone. 

MODEL DEVELOPMENT 


A three-dimensional model of the nearshore circulation is derived assuming 
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stationary wave conditions, straight and parallel bottom contours, and random waves that are 
narrow-banded in both frequency and direction. A right-handed-coordinate system is used 
with x positive onshore, y positive to the south and z positive upward from the sea surface. 
In the following, the energy flux balance is considered first to describe wave transformation, 
which includes contributions by breaking waves described with the roller concept. The 
random wave distributions are described, which are used to ensemble average the various 
equations. The conservation of mass is described next, which provides an integral condition 
for the solution of the vertical profile of the cross-shore velocity. The conservation of cross- 
shore and alongshore momentum fluxes are considered next. The depth integrated cross- 
shore momentum equation describes wave set-up/down, which is the primary forcing for 
cross-shore currents (undertow). The alongshore momentum equation contains the lateral 
transfer of mean momentum term that couples cross-shore and longshore currents, which 
requires specifying the vertical profile of U(x,z) by solving the cross-shore momentum 
equation. The depth integrated alongshore momentum equation is then solved to find the 
cross-shore distribution of the mean longshore current, V(x). 
Conservation of Energy Flux 

Wave heights across the surf zone are determined using the energy flux balance 


equation 
“ [EC cos(a)]} + ai [E_Ccos(a)] = D 1 
oc. Oe (1) 


where the energy 1s partitioned into wave, E_, and roller, E, contributions (Svendsen, 1984), 


C, and C are the wave group velocity and phase speed described by linear theory, « is the 
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mean incident wave angle with respect to the shore normal calculated using Snell's Law, and 
D represents dissipation of roller energy (Nairn er al., 1990). Following Deigaard (1993), 
the roller dissipation is modeled as D = C ts where o is the mean surface shear stress at 
the wave/roller interface. 

Both cross-shore distributions of wave height and the roller dissipation, D, are 
calculated using the Lippmann and Thornton (1997) model. This model assumes that random 
wave heights are Rayleigh distributed everywhere, even when breaking (Thornton and Guza, 
1983) and was shown to give good results for the DUCK94 data set (Garcez Faria et al., 
1997b [Chapter 2]). Therefore, wave properties in the various equations are ensemble 
averaged by integrating through the Rayleigh distribution, p(H), and using a weighting 
function, W(H), (Whitford, 1988) to calculate the percentage of waves that are breaking, 
DAH) = WH) ptf). 

Conservation of mass flux 


The conservation of mass for straight and parallel contours is given by 


n 
| - p[ U(z) + H(z) +u(z)]dz = 0 (2) 


=i 


where the boundary condition of no flow through the beach has been utilized and the 
velocities have been partitioned into mean, wave and turbulent contributions. In an Eulerian 
reference frame, there is an onshore mass transport by waves and an additional transport by 
the rollers of the breaking waves (Svendsen, 1984) limited to an upper region between the 


crest and trough. The ensemble and time-averaged onshore mass transport in the upper region 
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including contributions from wave rollers is balanced by a mean return flow below the trough 


(undertow) 


M1 
<q.> = -< i pU(z)dz> = -pU_h, (3) 
~h 


where 1, is trough elevation, U_ is the depth-averaged retum flow, h, is the depth below 
the trough, and <>indicates ensemble averaging (see Garcez Faria et al., 1997b [Chapter 
2], for details). 

Conservation of momentum flux 


The depth integrated and time averaged cross-shore momentum equation is given by 


ee! ainda (H+ Ayn (4) 
+ - + = = +h)— 
Ox Ox ou Ox cae Ox 











where 7 is the mean wave setup (down), Ss is the wave radiation stress associated with the 
momentum flux due to the wave motion, the second term on the lhs is the cross-shore 
gradient of the momentum flux associated with wave rollers, and the last term on the ths is 
the convective acceleration of the depth-averaged undertow. 


The depth integrated and time averaged alongshore momentum equation 1s given by 





aS A n, Ne - = 
ee UVdz + i+u)(V+v+vde )] = -t? + t” 
aes ie 50 a) ( dz J = -t + 73) (5) 


where F\, is the ensemble averaged alongshore wave forcing, which is modeled here using 
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a two-layer concept where waves are breaking (Lippmann er al., 1996). The lower layer 
contains the organized wave motion and the upper layer describes wave breaking as an 
elevated body of turbulent fluid riding on the front face of the wave (wave roller), consistent 
with spilling breakers. Following Deigaard (1993), the exchange of momentum across the 
wave/roller interface is described by the upward and downward volume fluxes between these 
two layers. The downward flux produces a shear stress at the surface, T., that 1s associated 
to the total wave forcing, F y? by conservation of momentum. Therefore, following Stive and 
De Vriend (1994), 1 can be written in terms of the average shear stress at wave/roller 


interface and local mean incident wave angle 


p « Soe , Sor. = sinh 
a wot = Cesc 
Ox Ox ; (6) 








where ee is the alongshore wave radiation stress and See, is the roller contribution to the 
radiation stress. The surface shear stress is calculated from the roller dissipation, D, using 
Lippmann and Thornton (1997) model, which is calibrated with video records of the cross- 
shore variation of wave breaking (see Garcez Faria et al., 1997b [Chapter 2], for details). 
The second term on the lhs of (5) is the lateral mixing associated with the integrated 
turbulent Reynolds stress that is parameterized using an eddy viscosity model (Bowen, 1969; 


Thornton, 1970) 


— — @V 
S.. = ak = Bee) (7) 


13] 


where v, is a kinematic eddy viscosity that is estimated following Battjes (1975) by 
Vea Cia n) 2)", where N is an adjustable coefficient of O(1), the turbulent length 
scale is order of the local depth, and the turbulent velocity scale is Ay" 

The third and fourth terms on the lhs of (5) represent the lateral mixing associated 
with the interaction between cross-shore and alongshore currents (Putrevu and Svendsen, 
1993), which require specifying the vertical profiles of both longshore and cross-shore 
currents and will be addressed later. 

Following Thornton and Guza (1986), the time and ensemble averaged non-linear 


alongshore bottom shear stress, tT), is modeled as 


<t’> - [tz feg (V?2+a?+2Vasin(a) (V+asin(@)) dt] p()dH — (g) 


where C f is a bed shear stress coefficient. 
The time averaged alongshore wind stress, T” , is calculated using the drag coefficient 


method 


tT’ = pC, WW, (9) 


where p, is the atmospheric density, W and We are the total and alongshore components of 
wind speed and C,, is an atmospheric drag coefficient. For C,,, the value obtained by the 
WANDI Group (1988) from observations at intermediate-depth water, 1s enhanced by 33 
percent to take into account the increased surface roughness associated with breaking waves 


(Whitford and Thornton, 1993, 1996). 


eo 


Vertical profile of longshore velocities 

Garcez Faria et al. (1997a [Chapter 1]}) using data from DUCK94 showed that the 
vertical structure of mean longshore currents on a barred beach is well described by a 
logarithmic profile 


V + 
Viz) = — ee) (10) 
K La 


a 


where « is the von Karman constant (0.4), v_is the alongshore shear stress velocity, and z, 
iS an apparent roughness due to the increase in effective roughness owing to the presence of 
waves (Grant and Madsen, 1979). The vertical profile of the longshore current is 
parameterized by v_ and z,, and the mean bottom shear stress, is related to the alongshore 
shear stress velocity by T, oe 


To compare model predicted longshore currents with measurements, the mean 


longshore current found by integrating (10) over depth is defined as the reference velocity 





given by 
vay ed pf vein Zt 4) de = ve piney + 22-1] 
x) =— aT 2 Nt == ) 
4 a a K Z, h (11) 


Vertical profile of cross-shore velocities 
Garcez Faria et al. (1997b [Chapter 2]) extended the monochromatic wave 
formulation for undertow by Stive and Wind (1986) using the random wave formulation of 


Thornton and Guza (1983). The vertical structure of the ensemble averaged undertow Is 
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obtained by depth-integrating the cross-shore momentum twice and solving for the 
integration constants by applying conservation of mass over the vertical (2) and the stress at 
the trough level (by integrating the vertical momentum equation from the bottom to crest 


elevation), to give 


2 3h? - 
=U(z)> —- Us = vas SEEN] frp + 
pu 2 








<te 
ee N, 
es] poses 
[z a (12) 


where pu is a depth-independent vertical eddy viscosity coefficient (Garcez Faria et al., 


1997b [ Chapter Piet is the mean cross-shore bed stress, and F’ is the undertow forcing, 





given by 
eu” leo 7 a én 
he + ——|[u“-w<] + pg— 13 
ets ae RS te 
The cross-shore bottom stress, has already been shown to be small (Longuet-Higgins, 


and Stewart, 1964), and thus the last term on the rhs of (12) is neglected here. 
Momentum transfer by mean currents 

Momentum mixing is the result of the vertical shear between the cross-shore and 
longshore mean currents as described by the third and fourth terms on the |hs of (5). In the 
water column beneath the trough, the alongshore momentum of the mean longshore current, 


o V, is advected offshore by the undertow, giving 


n, 
belie - pis le analy i) - + 





h; S10 (14) 
KH, 
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which can be rewritten as a function of the reference velocity, V_ (11), neglecting the term 


oF : 
—«1,to give 
h 


n, 
[ pUVdz = p<U >h,V_, + 
-h 


v 


hy <F,> (15) 
KH, 





In the surface region between the crest and trough of the wave, the alongshore 
momentum of the mean longshore current, p V, is advected shoreward due to the mean mass 
transport velocity of the waves and rollers of breaking waves, which is described by the 
fourth term on the lhs of (5). This term can be simplified assuming that waves and turbulence 
are uncorrelated, small incident wave angles and that the radiation stress contributions (p uv) 


in the crest-trough region are of higher order, giving 


ae 1, 
z Fe — Z x, ih 
fe +u)(V+v+v)dz ] = pV(n=0) [ @ +u)dz = <q,> od (16) 
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1, 1, 


The cross-shore distribution of the mean longshore current, V_(x), is calculated by 
solving (5) using the Newton-Raphson iterative method, applying a centered finite 
differencing scheme to rewntte the partial derivatives with the boundary condition of no flow 
through the beach. 

FIELD DATA 

Field data to test the model were acquired during the DUCK94 nearshore expenment 

at the U.S. Army Corps of Engineers Field Research Facility (FRF), Duck, North Carolina, 


in October 1994. For details of the experiment see Thornton et al. (1997), Elgar et al., 1997, 
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and Garcez Faria et al. (1997a (Chapter 1], and 1997b [Chapter 2]). The data presented here 
are from 10-12 Oct. when strong currents caused by a storm with predominant winds and 
waves from north were present and the vertical structure of both longshore and cross-shore 
currents have been studied by Garcez Faria et al. (1997a (Chapter 1], and 1997b [Chapter 
2]). During this period, observations spanned the entire surf zone and conditions approximate 
the assumptions of steady state and straight and parallel contours. 

Cross-shore transects of wave heights and currents were measured using a cross-shore 
array of two-component Marsh-McBimey electromagnetic current meters (ems, hereafter) 
and pressure sensors at each of 13 locations (See Elgar et al., 1997, for details). The vertical 
structure of the current was measured with a vertical stack of seven ems mounted on a 
mobile sled, which was deployed approximately 25 m to the south of the fixed cross-shore 
array. Additionally, directional wave spectra were acquired using a linear array of 10 pressure 
sensors in 8 m depth. 

For the first run on each day, the sled was towed by the Coastal Research Amphibious 
Buggy (CRAB) to its farthest offshore location seaward of the bar (approximately 160 m 
from the shoreline). A forklift on the beach pulled the sled shoreward 10 to 30 meters for 
subsequent measurement runs that are referred to in the text by sequential numbers within 
each day. Each data run was nominally one hour, and seven to eight runs were made across 
a transect during each day spanning the high tide during this period. 

For the cross-shore array, longshore currents were measured near mid-depth, but at 
arbitrary elevations, z,.. To make comparisons with the model, measured longshore 


currents, V 


meas? are corrected to correspond to the reference mean longshore current, V_, 
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using 
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where V__ is the corrected velocity. The elevations of the ems relative to the mean water 
depth, z,__.,. were adjusted using bathymetric profiles measured daily by the CRAB and the 
measured mean tidal elevation. The apparent roughness, z,, at each sensor location within 
the cross-shore transect occupied by the sled was obtained by interpolating the values 
obtained at the sled locations found by fitting a logarithmic profile to the measured vertical 
distribution of longshore currents based on a linear-regression least-squares method (Garcez 
Faria et al., 1997a [Chapter 1]). For sensors located off-shore of the sled transect, a constant 
value equal to the measured zat the most off-shore position occupied by the sled was used. 
An example of the measured and corrected mean longshore velocities during the first run of 
11 Oct. are shown in the upper panel of Fig. 2. The root mean square (rms , hereafter) value 
of the velocity corrections for all 22 runs considered is 7.1 cm s” which results on a rms 
reduction of 16 percent, and maximum correction of 16.4 cm s' amounting to a 27.5 percent 
reduction (Fig. 2, lower panel). Data analysis showed that the velocity correction 1s not a 
strong function of z,, but rather mostly determined by the ems’ elevation and local depth. 
If instead of an interpolated value, a constant z, equal to the mean of all measured values 
within the days being analyzed were used, the rms difference between the velocity 


corrections given by these two approaches is 2.1 cm s’' (5 percent). 


iy 


RESULTS AND DISCUSSION 
The relative importance of the various terms in the depth integrated and time 
averaged alongshore momentum balance (5) is investigated by comparing model predictions 
with field measurements. Contributions from each term in (5) are addressed separately, first 
by using a constant bed shear stress coefficient (C >” and then by specifying cross-shore 


variation of C f° Model/data comparison are quantified using a relative rms error defined by 


e = 100[ (18) 





In model comparisons with data, the vertical structure of mean currents 1s examined 
at the single sled location during each run, while the cross-shore distribution of mean 
longshore velocity is examined using the array data combined with the single sled data point. 
The runs used in the following figures to illustrate contributions from different terms in the 
alongshore momentum balance are arbitrarily chosen to cover as much of the data as possible 
Alongshore wave forcing 

The importance of including contributions from wave rollers is examined in a simple 
balance between wave radiation stress and bottom shear stress using the Lippmann and 
Thornton (1997) model that includes rollers. The basis for comparisons is the Thornton and 
Guza (1986) model that does not include rollers, which implies that dissipation of wave 
energy by wave breaking occurs locally; this model has good success in describing wave 
driven longshore currents on near planar beaches. A comparison of the two model 


predictions with observations of mean longshore currents for the third run of each day being 
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examined is shown in Fig. 3. Incorporating the roller contribution into the alongshore wave 
forcing significantly changes the cross-shore distribution of the mean longshore current, 
resulting in a better agreement with observations (Tables 1 and 2). 

The effects of including wave rollers is to displace the forcing of the longshore 
current shoreward. The magnitude of this onshore displacement is controlled by the steepness 
of the front face of the breaker and wave phase speed (Nairn et al., 1990; Deigaard, 1993; 
Thornton and Lippmann, 1997). For highly asymmetric waves, the advection distances are 
smaller and the longshore current maximum will be located closer to the sand bar crest than 
for more symmetric waves. Advection distances also increase with increasing wave phase 
speed, which is proportional to water depth in shallow water, and hence more current is 
forced in the trough at higher tidal elevations. This tidal modulation of the longshore current 
is consistent with observations during the DELILAH expenment (Thornton and Kim, 1993). 

This good agreement with data suggests that the alongshore wave forcing within the 
surf zone can be modeled using the average shear stress at the wave/roller interface, oe as 
long as model predicted shear stress is constrained by wave breaking observations. 

It can also be noticed in Fig. 3 that good agreement is obtained between the sled data 
and the depth-corrected data from the fixed array, despite the approximate 25 meters 
separation in the alongshore direction between the sensors. These results support the 
assumption that there was no alongshore variation of the mean longshore current in the 
vicinity of the measurements during the period being analyzed. 

Wind forcing 


During the three-day period studied here the mean alongshore wind speed was 9.1 
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ms! (range 7.7 - 12.6ms’) resulting in an average wind force to wave force ratio of 0.21 
(0.10 - 0.50) [Chaper 1]. Despite the wave force dominance characterizing wave-driven 
longshore currents, wind force contribution 1s not negligible. The effect of including the time 
averaged wind stress contribution, ie to the alongshore momentum balance is shown in the 
upper panel of Fig. 4 for the time period with largest wind force to wave force ratio (seventh 
run of 10 Oct.). Better agreement with observations offshore and within the trough is 
obtained, which is related to the relative increase in the wind force to wave force ratio as 
wave forcing is decreased over these regions. 

Ge is a direct function of the wind drag coefficient, C,,, that is yet poorly constrained 
by measurements within the surf zone. The constant value for this parameter adopted here 
is questionable, as the increased surface roughness over the bar associated with wave 
breaking is theoretically expected to result in larger values for C, over this region. The 
sensitivity of the model to uncertainties in C, is shown in the lower panel of Fig. 4 by 
varying this value within the range of values obtained by Hsu (1970) for the upper swash 
zone of a Florida beach(C, = 1.0-7.0x 10°). These results indicate that uncertainties in C d 
could lead to underestimation of the wind stress contribution to the alongshore momentum 
balance. Hence, under certain conditions tT" can be an O(1) term even for wave dominated 
conditions inside the surf zone, especially within the trough of a barred beach. 

Turbulent Reynolds stress mixing 

Little difference is found over most of the cross-shore transect covered by the 

measurements between model predictions with and without the inclusion of the turbulent 


Reynolds stress mixing, Sm parameterized using an eddy viscosity model(Fig. 5, upper 
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panel). The large differences close to the shore are not of a physical nature, but rather 
limitations of the dissipation model as the depth approaches zero at the shoreline (see 
Lippmann and Thornton, 1997 for details). 

The gradual change in the predicted energy dissipation given by the Lippmann and 
Thornton (1997) model (Fig. 5, lower panel) results in an already smooth longshore current 
distribution that does not require this ad hoc mixing to smooth current gradients. Another 
possible reason that this turbulent mixing term has little effect on the longshore current 
profile is the small turbulent mixing scale by breaking waves used, O(h) or O(H,__), which 
is much less than the length scale of the width of the surf zone during the days considered 
050-100 h). Other turbulent mixing scales order of the surf zone, such as shear instabilities, 
may contribute to momentum mixing (Dodd et al., 1997), but are beyond the scope of this 
paper. 

Following these results, and to allow an independent evaluation of the importance of 
the momentum mixing by mean currents to the alongshore momentum balance, the turbulent 
Reynolds stress mixing is neglected for the remainder of the paper, by setting N to zero. 
Momentum mixing by mean currents 

The contribution from the momentum mixing by mean currents to the alongshore 
momentum balance is illustrated by comparing model predictions including and neglecting 
this term with observations of mean longshore currents for the second run of each day being 
examined (Fig. 6). The inclusion of this type of mixing results in a shoreward displacement 
of the peak of the predicted mean longshore current improving the agreement with 


observations within the surf zone for either a constant or a variable C f by O[10 percent] 
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(Table 2). Outside the surf zone, where both the magnitude and vertical structure of the 
cross-shore current decrease, this mixing effect is not as important, although it still leads to 
significantly smaller errors in the case of a constant C f (Table 1). It can also be noticed in 
Fig. 6 and Table 2 that better agreement with data is obtained towards the end of the three- 
day-period. It is postulated that this improvement in model/data comparison is due to 
enhanced non linear mean current interactions associated with the observed strong cross- 
shore currents caused by increased wave heights and smaller angles of incidence with respect 
to beach normal connected with the arrival of a storm (Garcez Faria et al., 1997b [Chapter 
2)). 

The solution for the vertical structure of the undertow used in describing the 
momentum mixing has as a free parameter the vertical eddy viscosity coefficient, u , that was 
kept constant (u = 0.04) in both the cross-shore and over the vertical, which was shown to 
be a reasonable approximation under field conditions (Garcez Faria et al., 1997b [Chapter 
2]). The robustness of this approach is shown by changing the values of by an order of 
magnitude (Fig. 7). The amount of mixing is directly proportional to and the model is 
stable for increasing values of u. However, reducing u by a factor larger than 5 leads to 
numerical instability as the vertical profile of the undertow (12) is an inverse function of u. 
These results indicate that the dispersion mixing is enhanced by an increase in the vertical 
gradient of the undertow (decrease in pu), which is consistent with theoretical expectations, 
and that the model is numerically stable for a wide range of values of the free parameter nu. 
Bed shear stress 


The alongshore forcing by winds and waves, including the effects of wave-breaking 
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modeled with the roller concept, plus the momentum mixing associated with non linear mean 
current interactions is balanced by the bed shear stress, modeled with a non linear quadratic 
formulation (Thornton and Guza, 1986). Best fit values for the constant bed shear stress 
coefficient, C rf? listed in Table | are found by iteratively minimizing the error, ©, between 
measurements and the complete model predictions (waves, rollers and UV mixing). Only 
ems located shoreward of the mean breaker line (at approximately 150 m offshore) are 
included in the error estimate in order to avoid biasing the rms fit by large errors from ems 
outside the surf zone. Nevertheless, the values of elisted in Table 1 are calculated using all 
ems in order to compare the effects of incorporating a variable C f° Due to the non-linear 
formulation used in the model, the values of e for partial solutions of the model (waves and 
rollers) are obtained by re-running the model until the C value that gives a rms best fit is 
found for each case. The average values of efor the ems only within the surf zone are listed 
in Table 2. 

Although a constant Cr value provides a reasonably accurate and robust numerical 
solution for the cross-shore distribution of mean longshore currents within the surf zone, it 
does not provide insight into the physics of the problem, nor foe it describe the circulation 
outside of the surf zone well (Figs. 3-7). Feddersen et al. (1997), using measurements during 
the entire DUK94 experiment balanced cross-shore integrated wind and wave forcing with 
a quadratic bottom stress formulation and found that best fit C f values within the surf zone 
are three times larger than the values outside of the surf zone. Garcez Faria et al. (1997a 
[Chapter 1]) equated bed shear stress calculated from logarithmic velocity profiles to a 


quadratic bottom shear stress formulation and found that the associated C f values varied by 
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more than an order of magnitude across the surf zone (0.0006-0.012) and were positively 
correlated with 7ms bottom roughness estimated from sonar altimeter measurements. Their 
measurements show that smaller values of C f occurred over the bar associated with small 
bottom roughness owing to the planning action of the intense wave-breaking in this region. 
Larger values of C y were measured in the trough, where the roughness increased due to the 
presence of mega ripples aligned with the strong observed longshore currents. They found 
the following empirical relation between C f and apparent roughness length scale of the bed, 


= 30 'Za normalized by water depth. 


] 


k 
Cy = 10,000 ales (19) 


Here, a first attempt 1s made to estimate the effects of a cross-shore varying C on the 
spatial distribution of mean longshore currents by comparing model predictions using this 
simple relation with measurements in a diagnostic approach. In order to use (19), measured z_, 
values for the entire three-day-period were linearly regressed as a function of cross-shore 
position with small values over the bar and increasing shoreward. A constant z_equal to the 
furthest offshore regressed value is used for cross-shore positions offshore of the transect 
occupied by the sled, where no measurements were available (Fig. 8, upper panel). The 
points with large estimated errors (€ > 50 percent) in either z, or C , were not included in 
the regression. These regressed values are then used in (19) to estimate the average cross- 
shore variation of C f within the period being analyzed (Fig. 8, lower panel). Despite the large 


data scatter resulting in low correlation between measured and regressed z, (r = 0.53), a 
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statistically significant correlation at the 95 percent confidence level ( r = 0.96) between 
predicted and observed C;, is obtained. This is clearly not the case if an average value for C f 
(dashed line) is adopted. 

The simplistic assumption that the cross-shore distribution of the bottom roughness 
length scale (k,) remains unchanged during the span of these three days which have 
significant changes in wave and wind forcing is certainly not realistic. Therefore, it was 
expected that the magnitude of the estimated cross-shore distribution of C f would have to 
be adjusted for each run by multiplying the results given by (19) by a factor of O(1). 
Nevertheless, no change was required for all the runs on 10 and 11 Oct., which suggests that 
the cross-shore variation of C f is not overly sensitive to changes in k,, and is mostly 
controlled by bathymetrical changes associated with tidal variation. On 12 Oct., all the runs 
had to be increased by a constant factor of 1.7 to match observations. This is not surprising, 
as the storm waves on this day increased bottom roughness within the trough (Thornton et 
al., 1997) and enhanced the non-linear interactions between waves and currents, which are 
the two major factors defining the magnitude of k, (Grant and Madsen, 1979). Increased z, 
and C, values were measured on 12 Oct. (Fig. 8), although these anomalously large values 
were excluded from the regression due to their large uncertainties. 

A comparison between model predictions using the best fit constant C f and the cross- 
shore variable C f given by (19) is shown in Fig. 9 for the first run of each of the days 
examined. These results, together with smaller calculated e (Table 1) show that despite the 
crude approximation adopted here, a variable C f significantly improves the agreement with 


observations outside the surf zone, O[20 percent]. Even inside the surf zone, where the 
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constant C f value was adjusted to give the best rms fit to observations, a similar or even 
better agreement 1s obtained (Table 2). The overall improvement of model/data comparison 
obtained by applying the quasi-three-dimensional model with a cross-shore varying C f 
indicates that the bottom roughness length scale is an important factor controlling the bed 
shear stress. It also suggests that the empirical relation obtained by Garcez Faria et al. 
(1997a [Chapter 1]}), (19), provides a reasonable estimate for the cross-shore variation of C f 
under field conditions, although the prognostic use of this relation is still limited due to our 
present inability to predict the bottom roughness length scale. 
SUMMARY AND CONCLUSIONS 

The predicted cross-shore distribution of mean alongshore currents over a barred 
beach given by the quasi-three-dimensional model formulated here are compared with field 
observations acquired during the DUCK 94 experiment in an attempt to improve our physical 
understanding of the hydrodynamic processes governing longshore currents in the surf zone. 
The model includes forcing due to breaking waves described using the roller concept 
(Lippmann and Thornton, 1997), alongshore wind stress, cross-shore advection of mean 
momentum of the alongshore current (Svendsen and Putrevu, 1994), and a full non linear 
bottom shear stress with a variable bed shear stress coefficient, C ip constrained by 
observations (Garcez Faria et al., 1997a {Chapter 1}). 

Incorporating the roller contribution results in a shoreward displacement of the 
forcing of the longshore current controlled by the wave asymmetry and phase speed. Using 
the shear stress at the wave/roller interface, calibrated with observations of the spatial 


variation of wave breaking (Lippmann and Thornton, 1997), to force the longshore current 
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significantly shifts the longshore current maximum towards the trough of the barred beach, 
resulting in largest improvement with observations, O[50 percent], for the various 
mechanisms examined. 

Improved agreement with observations offshore and over the inner trough of the sand 
bar is obtained by including the time averaged wind stress contribution, which is modeled 
using the drag coefficient method. This better agreement is related to the relative increase in 
the wind force to wave force ratio as wave forcing decreases over these regions. 
Measurements of wind stress within the surf zone are needed to better constrain the 
magnitude of the wind drag coefficient,C,,, and to investigate the expected cross-shore 
variation of this parameter with increased surface roughness associated with breaking waves. 

Two mixing mechanisms for the mean longshore current were investigated. First, the 
turbulent Reynolds stress mixing 1s modeled using an eddy viscosity formulation (Battjes, 
1975), and it is shown that it does not significantly improve the agreement with observations. 
The failure of this mechanism is possibly related to the scale used to parameterize this 
turbulent mixing by breaking waves, O(h) or O( 7, .), which 1s much less than the length 
scale of the width of the surf zone during the days considered O(50-100 h). The second 
mechanism, the momentum mixing by non linear mean current interactions, 1s modeled by 
combining the solutions for the vertical profiles of longshore and cross-shore currents 
Rearce7 Faria et al., 1997a [Chapter 1] and 1997b [Chapter 2]). This mechanism is shown 
to displace the peak of the predicted longshore current shorewards resulting in a better 
agreement with observations within the surf zone O[10 percent]. 


The use of anon linear bed shear stress formulation with a cross-shore varying bed 
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shear stress coefficient, C constrained by observations significantly improves the agreement 
with observations, O[20 percent], compared with using a constant C f value adjusted to give 
the best rms fit to observations inside the surf zone. These results indicate that incorporating 
the effects of a changing bottom roughness length scale through a variable C f iS an Important 
step for understanding nearshore hydrodynamics. However, a predictive relation for the 
bottom roughness length scale as a function of mean currents in the presence of waves is 
required to accurately specify C f values. 

The overall good agreement obtained between model predictions and observations 
indicates that alongshore pressure gradients were weak during the period being analyzed (10- 
12 Oct.), and hence the circulation was essentially two-dimensional. Therefore, despite the 
lack of setup measurements, the DUCK 94 data provide a detailed data set to verify existing 
2D and quasi-3D models for longshore currents under field conditions. 
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Figure 1. Cross-shore distributions of predicted (Lippmann, Brookins, and Thornton, 1996 
model) H_ wave heights (upper panel) and mean longshore currents (Thornton and Guza, 
1986 model) V_ (lower panel). The open circles are the measured values on 10 October 1994 


during DUCK94 experiment. The bottom profile (not to scale) is shown in the upper panel. 
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Figure 2. Cross-shore distribution of measured (+) and depth-corrected (0)mean longshore 
currents, V_, for the first run of 11 Oct. (upper panel). Cross-shore distribution of the 
calculated velocity-correction using (17) applied to the measured longshore velocities in the 
fixed cross-shore array for 10-12 Oct (lower panel). The bottom profile (not to scale) is 


shown in the upper panel. 
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Figure 3. Comparison between measured mean longshore velocity at the sled (0) and fixed 
array (+) and predictions forcing the model with waves only (solid line) and waves and 
rollers (heavy line) for the third run of 10 Oct. (a), 11 Oct. (b), and 12 Oct. (c). Respective 


bottom profiles (not to scale) are shown in all panels. 
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Figure 4. Comparison between measured mean longshore velocity at the sled (0) and fixed 
array (+) and model predictions neglecting (solid line) and including (heavy line) 
contributions from the alongshore wind stress for the seventh run of 10 Oct. (upper panel). 
Model sensitivity to the wind drag coefficient, C, (lower panel). The bottom profile (not to 


scale) is shown in the upper panel. 
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DUCK94 — Reynolds Stress Mixing 
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Figure 5. Comparison between measured mean longshore velocity at the sled (0) and fixed 
array (+) and model predictions neglecting (solid line) and including (heavy line) the 
Reynolds stress mixing for the fifth run of 11 Oct. (upper panel). Cross-shore distribution of 
predicted dissipation given by the Lippmann and Thornton (1997) model (lower panel). The 


bottom profile (not to scale) is shown in the upper panel. 
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Figure 6. Comparison between measured mean longshore velocity at the sled (0) and fixed 
array (+) and model predictions neglecting (solid line) and including (heavy line) the 
momentum mixing by mean currents for the second run of 10 Oct. (a), 11 Oct. (b), and 12 


Oct. (c). Respective bottom profiles (not to scale) are shown in all panels. 
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Figure 7. Model sensitivity to the vertical eddy viscosity coefficient, u, compared with 
measured mean longshore velocity at the sled (0) and fixed array (+) for the first run of 11 


Oct. superposed on the bottom profile (not to scale). 
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Figure 8. Apparent roughness, z, , versus cross-shore distance (upper panel). The dashed line 
represents the regressed values used in (19). Bed shear stress coefficient, C fo Versus CIOSs- 
shore distance (lower panel). The solid line represents predictions given by (19) using 


regressed z, values shown in the upper panel and the dashed line is the average C f° 
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Figure 9. Comparison between measured mean longshore velocity at the sled (0) and fixed 
array (+) and model predictions using a constant bed shear stress coefficient, C f° (solid line) 
and a variable C (heavy line) tor the first run of 10 Oct. (a), 11 Oct. (b),.and 12 Oct_(c). 


Respective bottom profiles (not to scale) are shown in all panels. 
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Table 1. Best fit model parameters and percent errors for all current meters 


Constant C; Variable C; 
Day Run Percent Meron PErceni Eton 


UV Average UV 
Cr Waves Rollers Mixing Cr Waves Rollers Mixing 


10 1 0.0050 63 50 51 0.0031 54 43 36 
10 2 0.0065 65 69 50 0.0030 65 44 43 
10 3 0.0050 60 43 45 0.0029 95 37 37 
10 4 0.0065 70 46 45 0.0030 ~—- 65 38 3G 
10 5 0.0060 276 52 50 0.0030 68 47 46 
10 6 0.0060 72 50 50 0.0032 54 37 38 
10 7 00055 85 55 48 0.0032 51 36 38 
11 1 0.0045 98 62 50 0.0031 57 23 23 
11 2 0.0055 67 41 40 0.0030 52 30 33 
11 3 0.0050 68 44 44 0.0030 54 33 33 
11 4 0.0060 74 47 46 0.0030 61 37 ay 
11 5 00070 78 54 48 0.0030 62 43 44 
11 6 0.0065 90 61 45 0.0030 ~—-62 a2 33 
11 7 0.0060 65 58 54 0.0032 54 37 41 
11 8 0.0090 83 64 59 =—s-0.0031_— 63 35 Ue 
12 1 000900 74 °}#&53 41 0.0030 62 33 23 
12 2 00085 74 45 20 0.0053 = «53 30 28 
12 3 0.0090 66 54 35. 0.005256 37 29 
12 4 0.0065 64 48 35 0.0050 56 29 on 
12 5 0.0070 66 44 40 0.0050 60 36 36 
12 6 0.0090 72 62 46 0.0050 62 48 45 
12 7 0.0100 68 139 98 0.0052 67 137-105 

Mean 0.0068 = 82 58 47 0.0036—Ss«CO6 41 39 
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Table 2. Average percent errors for current meters located within the surf zone 


Constant C; Variable C; 
Day Percent Error Percent Error 


Waves Rollers UV Waves Rollers UV 


Mixing Mixing 
10 110 47 oe ve 42 38 
el 89 ay 34 65 40 34 
12 78 a ae 65 By 34 
Mean 92 50 33 68 40 55 


163 


164 


INITIAL DISTRIBUTION LIST 


Defense Technical Information Center 
8725 John J. Kingman Rd., 


¢* «© e& e e & e& & © © © &© © © © © © © © © © © © © © © © © © &© © & © 


STE 0944 


Ft. Belvoir, VA 22060-6218 


Dudley Knox Libary 
Naval Postgraduate School 
411 Dyer Rd. 

Monterey, CA 93943-5101 


Oceanography Department 
Code OC/Co 

Naval Postgraduate School 
833 Dyer Rd. Rm. 328 


Prof E. B. Thornton 
Oceanography Department 
Code OC/Tm 

Naval Postgraduate School 
833 Dyer Rd. Rm. 328 


>. ee 8 ¢ @ 


LCDR A. F. Garcez Fana 


ee @ e# @ @ # @ @ @ @ @® © @® ¢ @ © @® @® @ @® ® &® ® ® ® ® @ @® @© @® 8 @ © © @© © © &© &© @ 6 


eoeee ee & e# @ © © © © @ © © #& © © © © © © © © © & © &© © &© © © © © © © © © © © & © 


Ss ec c+ es 8 e+ eC ee ee ¢ 8 Se ee ee een en ee Se eS hUmSehUhmDhUcMLUlUlUCehlUlUSlUlUe he! 6S Lt os 8 6. eh Ue hCe 


Av. Governador Portela, 810 


Manejo 


Resende, Rio de Janeiro, 27520-250 


Brazil 


Diretoria de Hidrografia e Navegacao 


Rua Barao de Jaceguay s/n 
Ponta da Armacao 


ee e«eeee@ e@ @®8 e@ #&#® @ @ @ @ © @ © @ © ® @ @ @ ® @® @® @® @ @® @® @ @ @© @ 


Niter6i, Rio de Janeiro 24048-900 


Brazil 


Prof. T. Stanton 
Oceanography Department 
Code OC/St 

Naval Postgraduate School 
833 Dyer Rd. Rm. 328 


«ef @ © @ @ & @ 


Proi. 72 HeC. Herbers; gainer, 6 ee oe 
Oceanography Department 

Code OC/He 

Naval Postgraduate School 

833 Dyer Rd. Rm. 328 


Ope J ALIp pin ant ec ta ees oe renee oe es co, 
Center for Coastal Studies, 0209 

Scripps Institution of Oceanography 

9500 Gilman Dr. 

La Jolla, CA 92039-0209 


CEO AUTO ARNG lene canteen soos SOM arcercecapeannett ace oe Pia tea tage ec aie ie 
Naval Oceanography Command 
Stennis Space Center, MS 39529-5000 


COMMMANGINGRO ICED |. 5 2.2 os wc bet es wa ns ee 
Naval Oceanographic Office 
Stennis Space Center, MS 39529-5001 


SUPCMINUCIMEIG ae. sais dises s eee tec aee ene tae ee a «2 ee 
Naval Research Laboratory 

Marine Geosciences Div Code 7400 

Stennis Space Center, MS 39529-5000 


@ivjetorinaval Research ..a.601siee4 500.5... 00 Ce ck eee 
800 N. Quincy Street 
Arlington, VA 22217 


OiMmiCe OrNavaleixeSealch >. 4090s. .ae%..... . . weet eee 
Ocean Sciences Directorate (Code 1121 CS) 

Attn: Thomas Kinder 

800 N. Quincy Street 

Arlington, VA 22217 


DET Tea eae e555 Sith ick. and. da ack pt os oe sdore cone lece «cs oe et ea 
Scripps Institution of Oceanography 
La Jolla, CA 92037 


Dr. Dennis Wy Miter)... 606 6 +0205 Bee bees. ee 
Oceanography Department 

U.S. Naval Academy 

Annapolis, MD 21402 


166 


lee: 


18. 


ND: 


Zt. 


o\, 


icc (Ons es es Se dee ke ee ee 


U.S. Army Coastal Engineering Research Center 
Kingman Building 
Ft. Belvoir, VA 22060 


ID) oe 1 ee 2 2 ce ek Re ek cee ee ee baw Wes 


Waterways Experiment Station 
Corps of Engineers 

3909 Halls Ferry Rd 
Vicksburg, MS 39180-6199 


Dn Diane FOSter aa ee oe eos ore oe bees ered 6 Ok des Soi ga hed 


Department of Oceanography 
Dalhousie University 
Hallifax, NS 

Canada B3H 4J1 


WDtecle lami auimelt-= «0.0 ced ene aks vo ee ee wre ce ew oe eee. Oe 


Institut fuer Meereskunde 
Universitaet Hamburg 
Troplowitzstr. 7 

D - 22529 Hamburg 
Germany 


MP ePAGIRENIETS s..44 « o..4eads Doe wba ds b4Gek se ele eben. ben Hoe. ee 


Delft Hydraulics 
Rotterdamseweg 185 
PO Box 177 

2600 MH Delft 

The Netherlands 


167 











